NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


DISSERTATION 


TIDAL  DYNAMICS 
AND  MIXING 

OVER  STEEP  TOPOGRAPHY 

by 

Emanuel  Ferreira  Coelho 
June  1994 

Dissertation  Supervisor  Timothy  P.  Stanton 

Approved  for  public  release;  distribution  is  unlimited 


94  8  93  031 


UNCLASSIFIED 

MtOjMTr  CUkSWCATKM  OF  THM  PA« 


REPORT  DOCUMENTATION  PAGE 


UNCLASSIFIED 


N/AVAILABIU 

Approved  for  public  release; 
distribution  is  unlimited 


6c  ADDRESS  (Oty.  Smm,  am/ZIPCodai 

Monterey,  CA  93943-5000 


OtMppkeactm) 

35 


Naval  Postgraduate  School 


7b.  ADDRESS  (Oty.  smm.  a/m  HP  Coat) 

Monterey,  CA  93943-5000 


■  7  ^  ^  •T ;!  :1IT71  *11  iT^T  *1  *  TT  J  fr 


1 1  i  a  [•?  *  I  -  I 


ORGANIZATION 

(ifatmteaua, 

8c.  ADDRESS  (City,  Smm.  ana  HP  Caaat 

ELEMENT  NO.  I  NO. 


NO. 

NO. 

1 1.  TTITjE  (mmutm  Sacunm  cmamtrrninn) 

TIDAL  DYNAMICS  AND  MIXING  OVER  STEEP  TOPOGRAPHY 


COSATI  COOES 


GROUP  I  SUB-GROUP 


18.  SUBJECT  TERMS  (OontM  an  /mymm  tfnaeaaamy  ant  manUy  by  Uoeknmbar) 

Internal  Tide,  Tldai  Mixing,  Internal  Waves,  Coastal  Dynamics,  Nonlinear 
Internal  Wave  Propagation,  Shelf  Break  Tidal  Dynamics 


19.  ABSTRACT  (Contnua  on  rmmaa  H/meaaamymttmanmybymocttnmba/i 

Coastal  currents  and  water  colunn  stratification  frequently  are  detexeined 
by  local  phenomena  that  interact  to  produce  complex  dynamics.  A  process  study 
was  designed  to  assess  the  role  of  tides  over  a  coastal  shelf  break  as  a  local 
source  of  internal  waves  spatially  inhomogeneous,  baroclinic  residual  currents, 
and  internal  mixing.  The  teiqioral  evolution  of  three  dimensional  stratification 
and  currents  were  characterized  within  a  4Kmii4Km  box  close  to  the  mouth  of 
Nonterey  Camyon,  Calif oznia.  Strong,  baroclinic,  unidirectional,  spatially 
inhoBogeneous  jet-like  currents  were  Observed  with  24  hour  periods,  together 
with  the  occurrence  of  semidiurnal  thermocllne  oscillations  and  vertically  mixed 
layers. 

The  linear  theory  of  Baines  (1982)  predicts  the  existence  of  significant 
tidally  forced  internal  waves  in  the  area,  with  an  Interfacial  mode  at  the 
thermocllne  decoupled  from  the  internal  waves  propagating  in  the  weaker 
stratification  below.  A  ray  tracing  method  is  developed  to  study  the  nonlinear 
'propagation  of  these  tidally  forced  internal  waves  by  considering  wave  advection 
and  straining  by  the  topographically  controlled  barotropic  tidal  flow  while 


tamon 


OOPORM  1473,84  MM 


83  APR  tdWen  may  b8  uaM  unS  •muBMd 
AH  otiar  adSona  am  obMW 


SECURTtV  CLASSIBCATION  OP  TH8  PAGE 

UNCLASSIFIED 


UNCLASSIFIED 

SECUWTV  ClASaif  ICATION  Of  THa  PAOC 

n«9l«cting  internal  wava-wav«  int«ractions.  It  is  found  that  nonlinear 
interaction  of  tidally  forced  internal  waves  with  the  topographically  controlled 
barotropic  tidal  flow  is  a  eechanism  capable  of  forcing  residual  currents  and 
internal  wave  shear  instabilities  that  are  capable  of  creating  internal  mixed 
layers  similar  to  those  observed  at  the  mouth  of  Monterey  Canyon. 


ii 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


Approved  for  pub 


release;  distribution  is  unlimited 


Tiual  Dynamics  and  Mixing 
Over  Steep  Topography 


by 


Emanuel  Ferreira  Coelho 
Lieutenant  Commander,  Portuguese  Navy 
B.S.,  Portuguese  Naval  Academy,  1983 
M.S.,  Naval  Postgraduate  School,  1991 


Submitted  in  partia' 
requirements  f 


ulfillment  of  the 
h-  degree  of 


DOCTOR  OF  PHILOSOPHY  IN  PH  IlCAL  OCEANOGRAPHY 

from  the 


NAVAL  POSTGRADUATE  SCHOCL 
June  1994 


Author: 


CRA&I  “Jj 

OTIC  TAB  A 

Unannounced 
JListification 


□ 


By . 

Dist.  ibution  / 


Availability  Codes 


Dist 


Avail  anu/or 
Special 


Approved  by: 


Emanuel  Ferreira 


Timothy  P.  Stanton 
Associate  Research  Professor 
of  Oceanography 


Iho 


James  H.  Miller 
jociate  Professor  of 
Electrical  and  Computer 
Engineering 


4: 


Roland  W.  Garwood 
Professor  of  Oceanography 


Robert  L. 
Professor  and  Chaii 
Meteorology 


lan  of 


Edward  B.  Thornton 
Professor  of  Oceanography 
Committee  Chalrmain 


ABSTRACT 


Coastal  currents  and  water  column  stratification  frequently 
are  determined  by  local  phenomena  that  Interact  to  produce 
complex  dynamics.  A  process  study  was  designed  to  assess  the  role 
of  tides  over  a  coastal  shelf  break  as  a  local  source  of  internal 
waves  spatially- inhomogeneous,  baroclinic  residual  currents,  and 
internal  mixing.  The  temporal  evolution  of  three  dimensional 
stratification  and  currents  were  characterized  within  a  4IQnx4IQn 
box  close  to  the  mouth  of  Monterey  Canyon,  California.  Strong, 
baroclinic,  unidirectional,  spatially  inhomogeneous  jet-like 
currents  were  observed  with  24  hour  periods,  together  with  the 
occurrence  of  semidiurnal  thermocline  oscillations  and  vertically 
mixed  layers. 

The  linear  theory  of  Baines  (1982)  predicts  the  existence 
of  significant  tidally  forced  internal  waves  in  the  area,  with 
an  interfacial  mode  at  the  thermocline  decoupled  from  the 
internal  waves  propagating  in  the  weaker  stratification  below. 
A  ray  tracing  method  is  developed  to  study  the  nonlinear 
propagation  of  these  tidally  forced  internal  waves  by  considering 
wave  advection  and  straining  by  the  topographically  controlled 
barotropic  tidal  flow  while  neglecting  internal  wave-wave 
Interactions.  It  is  found  that  nonlinear  interaction  of  tidally 
forced  internal  waves  with  the  topographically  controlled 
barotropic  tidal  flow  is  a  mechanism  capable  of  forcing  residual 
currents  and  internal  wave  shear  instabilities  that  are  capable 
of  creating  internal  mixed  layers  similar  to  those  observed  at 
the  mouth  of  Monterey  Canyon. 
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I.  IHTRODUCTIOM 


A.  MOTZVATZOM 

Coastal  currents  and  water  column  stratification  are 
determined  by  local  surface  forcing,  enhanced  tides,  rivers, 
topography  and  the  coastal  boundary,  together  with  larger 
scale  ocean  dynamics.  These  characteristics  determine  and 
control  circulation  patterns  and  distinguish  coastal  regions 
from  the  deeper  ocean.  Recently,  the  development  of  industries 
related  to  the  exploration  of  sea  resources  have  increased 
significantly,  and  methods  are  now  required  to  balance 
environmental  Impacts  and  industry  needs.  The  increase  in  the 
strategic  importance  of  these  areas,  where  most  of  the 
interaction  between  man  and  the  ocean  occurs,  should  introduce 
some  priority  to  understanding  the  dynamics  of  these  coastal 
regions . 

The  study  of  specific  mechanisms  contributing  to  advection 
and  mixing  in  coastal  areas  over  irregular  topography 
constitutes  the  broad  objective  of  this  research.  This 
objective  is  fundamental  to  a  large  number  of  fields. 
Including  the  study  of  sediment  transport,  pollutant 
dispersion,  coastal  air-sea  interaction,  marine  biology,  oil 
drilling,  coastal  engineering,  and  ocean  acoustics. 
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A  process  study  was  designed  to  assess  the  role  of  tides 
and  wind  forcing  on  mixing  and  advection  close  to  the 
continental  shelf  break,  where  steep  topography  and 
stratification  are  important.  A  3  day  surface  mixing 
experiment  was  executed  in  September  1991  to  assess  some  of 
the  complex  dynamics  in  Monterey  Bay.  An  Acoustic  Doppler 
Current  Profiler  (ADCP)  and  a  Tow-yo  CTD  package  was  used  to 
measure  currents  and  stratification  with  approximately  200m 
horizontal  resolution  along  a  closed  square  path  with  4Km 
legs,  sampled  approximately  every  2  hours.  The  measurement 
site  was  on  the  southern  side  of  the  Monterey  submarine 
canyon,  close  to  the  canyon  mouth,  over  highly  irregular 
topography  ranging  from  600  to  1500m  depth  (see  Figure  II- 
A.l) .  The  complex  dynamic  interactions  summarized  in  Figure  I- 
A.l  were  considered  during  this  study. 

A  brief  sximmary  of  the  experimental  observations  which 
motivated  this  work  follows.  The  winds  were  moderate 
throughout  the  obsexrvation  period,  resulting  in  a  very  shallow 
surface  mixed  layer.  Strong  spatial  inhomogeneities  in  the 
upper  ocean  current  and  temperature  stiructure  were  observed  in 
the  surveyed  area.  The  analyzed  oceanographic  data  suggest 
that  the  measurement  box  intersects  an  anticyclonic  mean 
baroclinic  flow  of  the  order  of  0.05ms'\  Strong  baroclinic  and 
spatially  inhomogeneous  northward  jet-like  features  were 
observed  24  hours  apart  in  the  south-west  portion  of  the 
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surveyed  area.  The  maximum  currents  occurred  between  50  to 
150m,  with  0.30ms'^  peak  velocities.  These  jet-like  pulses 
shallowed  with  time,  reaching  the  surface  at  the  end  of  the 
measurements . 

Thermocline  oscillations  at  tidal  periods  having 
amplitudes  of  up  to  20m  were  observed  at  the  north-west  corner 
of  the  surveyed  area.  The  internal  wave  field  was 
characterized  by  an  interfacial  mode  at  tidal  frequency  and  by 
smaller  scale  internal  waves  propagating  in  the  stratified 
layer  below  the  thermocline.  Vertical  thermocline  oscillations 
were  observed  above  the  current  maxima,  with  episodic 
vertically  homogenized  patches. 

An  internal  mixed  layer  was  observed  within  the 
thermocline,  between  20m  and  50m  depth  at  the  west  part  of  the 
surveyed  area,  with  strong  entrainment  of  warmer  water  down 
into  the  internal  mixed  layer,  occurring  in  an  area  of 
intensified  vertical  shear. 

In  the  present  work  it  will  be  shown  that  the  internal 
tides  (Baines,  1982;  Huthnance,  1989)  were  responsible  for  the 
observed  thermocline  oscillations  and  for  the  occurrence  of 
strong  vertical  shear.  Furthermore,  the  nonlinear  interaction 
of  the  internal  tides  and  the  topographically  enhanced 
barotropic  tide  (Willmot,  1987;  Maze,  1987;  Pichon  et  al., 
1990)  will  be  studied  using  a  numerical  model  to  show  that  it 
is  the  probable  cause  of  the  strong  baroclinic,  spatially 
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inhomogeneous,  unidirectional  currents  observed  at  the  mouth 
of  the  canyon.  Finally,  the  role  of  the  tidally  forced 
baroclinic  velocity  field  in  forcing  strong  mixing  events 
(New,  1988) ,  like  the  event  observed  during  the  MBAY 
experiment,  will  be  discussed.  This  work  is  organized  as 
follows.  Chapter  I  describes  the  motivation  and  the  scientific 
objectives.  Chapter  II,  describes  the  Monterey  Bay  Mixing 
Experiment  (MBAY)  conducted  in  Monterey  Bay  in  early  September 
1991,  and  the  data  processing  methods  used  in  the  analysis. 
Chapter  III  discusses  observations  motivate  the  hypotheses 
for  the  interpi  ■»tation  of  the  observed  phenomena.  Chapter  IV 
presents  the  theoretical  approaches;  the  linear  theory  of 
Baines  (1982)  is  used  first  to  estimate  the  magnitude  of  the 
forced  fields,  and  then  a  ray  tracing  method  is  developed  to 
study  the  propagation  of  the  tidally  generated  internal  waves; 
finally,  a  Richardson  number  criteria  is  used  to  estimate 
potentially  unstable  areas  where  mixing  is  likely  to  occur. 
Chapter  V  presents  the  conclusions  and  recommendations  for 
future  work,  emphasizing  the  scientific  contributions  of  the 
present  development. 

B.  SCIENTIFIC  OBJECTIVE 

The  oceanic  tide  is  a  barotropic  phenomena  that  can  be 
predicted  locally  (Gill,  1982).  However,  interaction  with  the 
coastal  topography  and  stratification  can  result  in  mean  tidal 
residual  flows  driven  by  potential  vorticity  conservation 
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Mchanisms  (Robinson,  1983;  Mass  et  al.,  1989a, b) ,  which 
contribute  significantly  to  changing  the  local  advection 
patterns.  Furthermore,  the  barotropic  tidal  motions 
interacting  with  the  upper  ocean  stratification  may  also  cause 
oscillations  in  density  gradients,  forcing  internal  waves 
groups  which  propagate  away  from  the  generation  areas  i . e . , 
internal  tides  (Baines,  1982).  These  tidal  forced  internal 
waves  are  often  non-linear  (Willmott,  1987;  Maze,  1987; 
Huthnanace,  1989) ,  and  may  force  mixing  events  by  vertical 
shear  instability  mechanisms  when  the  gradient  Richardson 
number  becomes  sufficiently  small  (New,  1988).  However,  none 
of  the  existing  theories  account  for  the  existence  of  the 
strong,  baroclinic,  unidirectional  current  pulses  obsezved 
during  the  MBAY  experiment. 

To  evaluate  the  importance  of  internal  tides  in  the  survey 
region,  the  tidally  forced  near-inertial  internal  wave  field 
is  estimated  using  the  linear  model  of  Baines  (1982).  The 
model  uses  idealized  topography  which  approximates  the 
continental  slope  close  to  the  southern  part  of  the  Monterey 
Bay,  and  is  forced  by  an  estimated  barotropic  tidal  field.  The 
modeled  velocity  field  is  then  compared  with  the  MBAY 
observations . 

To  study  near-inertial  internal  wave  propagation  on  a  non¬ 
stationary,  inhomogeneous  mean  flow,  a  nonlinear  numerical  ray 
tracing  method  is  developed,  neglecting  wave-wave  interactions 
but  retaining  the  advective  and  straining  effects  of  the 
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underlying  flow.  This  model  also  uses  idealized  topography 
approximating  the  continental  slope  in  the  southern  part  of 
Monterey  Bay,  and  is  forced  with  the  Baines  model  internal 
wave  field  at  the  shelf  break  and  by  an  estimated 
topographically  controlled  barotropic  tidal  flow.  The 
interpretation  of  the  numerical  results  emphasizes  the 
existence  of  rectified  baroclinic  flow,  which  is  comparable  to 
the  observed  strong,  baroclinic,  unidirectional  currents.  The 
modeled  velocity  fields  and  observed  mean  stratification  are 
also  analyzed  for  the  existence  of  potentially  unstable 
internal  shear  layers,  which  could  cause  vertical  shear 
instabilities  and  subsequent  mixing. 

The  present  work  is  based  on  the  internal  tide  propagation 
linear  models  of  Baines  (1982)  and  New  (1988)  ,  and  on  the  non¬ 
linear  studies  of  internal  tides  using  layer  models  by  Maze 
(1987)  and  Willmott  (1987),  where  only  tho  advection  was 
considered.  The  primary  contribution  of  this  work  is  an 
extension  of  Kunze  (1985),  who  studied  the  propagation  of 
near-inertial  internal  waves  on  geostrophic  shear,  by 
considering  the  time  variability  of  the  underlying  velocity 
field  through  which  the  waves  propagate. 


II.  MOMTERBY  BAY  MIZI1I6  BZPBRIIIENT 


A.  DATA  COLLECTION 

To  obtain  physical  data  and  to  study  the  photo  adaptation 
of  biological  activity  under  the  influence  of  advection  and 
diffusion,  the  Monterey  Bay  Mixing  Experiment  (MBAY)  was 
executed  between  4-7  September  1991  in  Monterey  Bay  under  the 
sponsorship  of  the  US  Office  of  Naval  Research,  the  Commander 
of  Naval  Oceanographic  Command  (CNOC) ,  and  the  Monterey  Bay 
Aquarium  Research  Institute  (MBARI) .  Two  ships,  the  RV  "Point 
Sur"  (sponsored  by  CNOC)  and  the  RV  "Point  Lobos”,  and  a 
Remote  Operated  Vehicle  (ROV)  from  the  MBARI,  were  used  as 
platforms  to  deploy  a  Tow-Yo  Conductivity-Temperature-Depth 
(CTD)  profiler  (5-80m  depth,  0.25m  bins),  a  CTD  profiler  (0- 
120m)  with  water  sampling,  two  Acoustic  Doppler  Current 
Profilers  (ADCP)  (7-263m  depth,  4m  bins) ,  a  laser  scanner  with 
a  high  resolution  video  camera  for  direct  bulk  diffusivity 
measurements,  a  Coherent-Doppler-Velocity  (CDV)  package  for 
microstructure  measurements  of  shear,  velocity,  pressure, 
temperature  and  salinity,  and  two  meteorological  packages  for 
the  measurement  of  the  atmospheric  conditions  and  bulk  air- 
ocean  fluxes.  Both  ships  used  Global  Positioning  System  (GPS) 
navigation  with  50m  accuracy.  Simultaneous  satellite  Advanced 
Very  High  Resolution  Radiometer  (AVHRR)  imagery  and  Coastal 
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Ocean  Dynamics  Application  Radar  (CODAR)  maps  were  also 
analyzed  to  determine  background  conditions. 

The  RV  "Pt.SUR"  followed  a  closed  square  path  with  4I0n 
legs,  repeated  27  times,  over  the  Monterey  Canyon  above  a 
highly  variable  ocean  floor.  ADCP  velocity,  temperature  and 
conductivity  profiles  and  meteorological  data  were  recorded 
with  spatial/temporal  resolutions  of  120m/60s  (ADCP  velocity) , 
360m/180s  (temperature  and  conductivity)  and  60m/30s 
(meteorological  variables).  The  RV  "Pt.Lobos**  was  stationed 
Inside  the  closed  path  recording  water  velocity  profiles  and 
meteorological  data,  and  controlling  the  MBARI  ROV  which 
measured  bulk  diffusivity  and  microstructure  levels  at 
selected  depths  between  10  and  50m.  These  measurements 
provided  direct  diffusion  and  turbulence  estimates  while 
simultaneously  characterizing  the  internal  wave  field,  tides, 
wind  stress,  and  surface  buoyancy  fluxes.  The  full  data  set 
provided  observations  of  mixing  in  surface  and  internal 
layers,  while  the  evolution  of  the  forcing  mechanisms  in  a 
sub-mesoscale  coastal  region  were  identified.  A  sketch  of  the 
survey  is  shown  in  Figure  II-A.l. 

Section  plots  were  made  of  the  ADCP  and  Tow-Yo  data  from 
the  RV  "Pt.Sur"  to  show  the  evolution  of  the  velocity 
components,  temperature,  salinity  and  water  density  profiles 
as  the  "Pt.Sur"  traveled  around  the  square  box  (one  plot  per 
box) .  Strong  jet-like  current  features  associated  with  layers 
of  strong  shear  occurred  periodically  throughout  the 
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observation  period  close  to  specific  depths  (surface,  upper 
themocline  and  lower  themocline) ,  as  shown  in  the  example  in 
Figure  II~A.2.  However,  these  strong  velocity  and  shear  events 
were  only  seen  at  the  same  locations  within  the  small  survey 
box,  demonstrating  the  strongly  inhomogeneous  character  of  the 
velocity  field. 

Considering  the  temperature  profiles  of  Figure  II-A.3,  the 
occurrence  of  isotherm  deepening  between  profiles  472  and  475 
within  the  strong  vertical  shear  region  (Figure  II-A.2) 
suggests  that  a  strong  mixing  event  occurred.  These  mixed 
layers  occur  well  within  the  pycnocline,  suggesting  they  are 
not  directly  forced  by  atmospheric  fluxes.  Understanding  how 
these  strong  shear  layers  appear  and  evolve  is  an  important 
step  in  characterizing  forcing  mechanisms  for  mixing  and 
diffusion  in  the  area. 

Well  resolved  temperature  inversions  were  frequently 
observed  within  the  thermocline;  however,  the  water  column  was 
still  hydrostatically  stable  as  salinity  compensation  is 
occurring  with  the  formation  of  cold  and  fresher  water  layers 
above  warm  and  saltier  water  layers  below,  as  the  T-S  diagrams 
of  Figure  II-A.4  show.  These  diagrams  are  displayed  for  two 
boxes  12  hours  apart  showing  a  general  reduction  in  density  in 
the  water  column,  particularly  at  30m  depth  where  sigma-theta 
decreases  from  25.7  to  25.2  with  the  water  becoming  warmer  and 
fresher.  The  generation  of  new  water  types  on  these  time 
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scales  can  only  be  explained  by  the  occurrence  of  strong  local 
nixing  or  by  lateral  intrusions  at  these  depths. 

The  COV  and  laser-scanner-video  data  from  the  ROV  were 
processed  separately  because  of  their  specific  characteristics 
and  scientific  purpose.  The  results  of  the  analysis  of  this 
data  will  be  reported  elsewhere. 

B.  DATA  PR0CB88ZM0 

After  the  preliminary  processing  and  analysis,  the  RV 
"Pt.Sur"  data  discussed  in  the  previous  section  were  spooled 
to  both  the  DOS  and  Unix  environments  for  further  processing 
using  MATLAB  software.  To  characterize  the  features  mentioned 
in  the  previous  section,  the  following  methodology  was 
employed: 

1.  Develop  an  Objective  Analysis  procedure  to  fit  the 
survey  data  into  a  regular  grid  to  facilitate  comparison 
between  different  data  types,  allow  uniform  averaging  of 
properties,  and  to  assess  the  spatial  and  temporal 
characteristics  of  the  observed  features. 

2.  Process  ADCP  data  using  a  reference  layer  to 
determine  relative  velocities  and  shear,  and  navigate  the 
reference  layer  using  the  ship's  Global  Position  System  (GPS) 
data  to  estimate  absolute  current  profiles. 

3.  Develop  a  spectral  estimator  to  detect  dominant 
propagating  waves  in  the  data,  based  on  the  repeated  box 

I 

J 
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sanpling  nethod  (time  series  following  similar  spatially 
closed  paths) . 


4.  Characterize  and  identify  strong  shear  layers  in 
space  and  time  from  the  ADCP  velocities  and  relate  them  to  the 
stratification  by  estimating  the  gradient  Richardson  numbers. 

5.  Characterize  the  internal  wave  field,  emphasizing 
near- inertial  frequencies  and  tide  components. 

6.  Identify  mixing  events  and  their  temporal/ spatial 
distribution  using  a  Richardson  number  criteria. 

7.  Perform  mass,  momentum  and  vorticity  balances. 

1.  Objective  Analysis  (OA) 

To  fit  the  ADCP  and  Tow-Yo  data  into  a  regular  time 
and  space  grid,  multi-dimensional  interpolation  methods  are 
used.  Ideally,  these  methods  are  based  on  correlation 
functions  obtained  from  the  data,  requiring  robust  statistics 
and  stationarity  not  seen  in  our  observations.  As  the  ADCP 
vertical  sampling  was  made  every  4m  and  the  Tow-Yo  sampling 
every  0.25m,  the  Tow-Yo  data  were  smoothed  and  then  decimated 
so  that  both  data  types  were  binned  every  4m.  Therefore,  the 
OA  was  performed  only  in  time  and  in  the  horizontal 
dimensions . 

Based  on  the  four  corners,  whose  coordinates  are  given 
in  Table  II-A-1.1,  four  legs  of  a  square  were  defined  with  120 
meter  spatial  sampling  and  1  minute  lags  between  consecutive 
samples  along  each  leg.  Several  analysis  methods  were  tested: 
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optiaun  interpolation,  linear  regression,  least  squares 
polynomial  fit  and  inverse  distance  weights.  Because  of  its 
form,  well-based  statistical  assumptions,  and  easy 
generalization,  the  optimum  interpolation  method  (Haltiner  and 
Williams,  1980;  Carter,  et  al.,  1986)  was  chosen.  In  this 
method,  it  is  assumed  that  the  field  at  the  grid  points  can  be 
obtained  from  a  linear  filter  operation  applied  to  the 
neighboring  observations  of  the  type: 

N 

^gxid  *  ^giid  “ 

where  the  Y's  are  the  field  values  (observed  and  at  grid 
points) ,  e's  are  the  error  and  the  p's  are  the  filter  weights. 
Following  a  least-mean-sguare  error  approach  and  assuming  that 
the  data  spans  a  stochastic  vectorial  subspace  where  the 
observed  inner  product  is  defined  by  the  operator  E [ ]  ( i . e . , 
ensemble  expected  value)  ;  the  least  square  condition  requires 
the  assimilation  error  vector  to  be  orthogonal  to  the  observed 
data  vectors,  i.e.,  error  and  measurements  exist  in  orthogonal 
subspaces.  Then 

-0  «> 

Defining  the  correlation  function  (CF)  as 

(3) 
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and  aasuaing  a  stationaxY  and  honogeneous  field,  the  function 
CF  can  be  nornalized  by  the  variance  of  Y.  The  solution  of  the 
least>square  condition  is  the  solution  of  a  hemitian  system 
of  equations,  with  ones  in  the  main  diagonal  and  smaller 
absolute  values  elsewhere,  in  the  form: 

(where  square  brackets  stand  for  a  matrix  and  underline  for  a 
vector) .  A  potential  problem  for  repeated  samples,  or  for 
purely  periodic  fields  is  that  they  may  generate  a  singular 
correlation  matrix.  A  solution  to  this  problem  is  to  identify 
repeated  or  periodic  samples  and  remove  them  from  the 
observations  before  the  analysis.  Another  approach  is  to  use 
an  eigenvalue  threshold  criteria  that  allows  handling 
redundant  observations,  or  ill  conditioned  matrices.  This 
procedure,  based  on  the  singular  value  decomposition  (SVD) 
(Strang,  1988;  Haykin,  1989)  allows  any  correlation  matrix  to 
be  inverted  and  gives  stable  numerical  results.  Defining  the 
pseudo-inverse  matrix  (PINV)  (Strang,  1988) 

PINV  »  0^^ Q  (5) 

where  [CF]  s  Q  E  Q^,  Q  is  the  orthonormal  matrix  of  the 
eigenvectors  of  [CF]  and  is  a  diagonal  matrix  with  the 
inverse  of  the  dominant  square  eigenvalues  of  [CF]  and  zeros 
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for  eigenvalues  smaller  than  a  threshold.  The  filter  impulse 
response  is  given  by 


£  =  [PINV\CE^  J  (6) 

and  the  mean  square  error  (MSE) ,  normalized  by  the  field 
variance,  is 

MSE  «  1-Ep*ce;  *  (7) 


As  the  limited  data  sampling  prevents  an  accurate 
estimate  of  the  two  dimensional,  depth  dependent  correlation 
function  of  each  variable  from  being  estimated  accurately,  the 
common  practice  of  using  an  appropriately  behaved  analytical 
correlation  function  has  been  implemented.  Based  on  the 
correlation  decrease  with  time  and  distance  for  non-periodic 
data,  several  functions  can  be  considered  and  tuned  for  each 
specific  application.  In  this  work,  the  function  displayed  in 
Figure  II-B-l.l,  similar  to  that  used  by  Carter,  (1986),  was 
used  with  450m  spatial  and  15  minute  temporal  scales.  These 
scales  are  based  on  the  physical  sampling  scales  and  are  not 
derived  from  the  observations.  This  function  was  used  in  the 
optimal  filter  to  weight  the  contributions  of  values  in 
surrounding  grid  points.  Valid  data  points  were  those  with 
distances  smaller  than  500m  and  with  time  lags  smaller  than  15 
minutes  from  each  grid  point.  Grid  points  with  no  valid 
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contributions  were  narked  as  bad  data.  The  RMS  errors  obtained 
from  this  filter  operation  are  an  indicator  of  the  tine/space 
distances  of  the  data  used  at  each  grid  point. 

The  "Pt.Sur"  navigation  was  irregular  with  frequent 
stops  in  order  to  obtain  CTD  casts,  and  sometimes  was  more 
than  500m  off  course.  This  fact  introduced  frequent  gaps  and 
additional  averaging  and  interpolation  was  required  close  to 
the  CTD  casts.  Maps  showing  the  vertical  averages  of  RMS 
errors  and  the  niimber  of  profiles  used  to  evaluate  each  data 
grid  point  for  each  box  are  displayed  in  Figure  II-B-1.2.  The 
frequent  stops  of  the  ship  appears  as  a  high  number  of  data 
points  at  a  given  location.  After  some  of  these  stops  the  RMS 
error  goes  to  one,  corresponding  to  a  data  gap  due  to  the  time 
threshold.  Nevertheless,  these  irregularities  in  the  seunpling 
rates  are  evenly  distributed  throughout  all  the  surveyed  area, 
so  no  significant  bias  is  expected.  However,  the  large  periods 
and  areas  with  no  data  can  introduce  some  local  aliasing, 
enhancing  or  attenuating  spatial  features.  Care  has  been  taken 
in  the  subsequent  analyses  to  avoid  using  boxes  with  large 
data  gaps. 

After  this  data  analysis,  linear  time  interpolation 
was  performed  at  each  grid  point  so  that  instantaneous 
estimates  of  data  fields  over  the  full  grid  were  available  for 
analysis. 
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2.  ADCP  V«looiti«a 

The  64  bln  ADCP  velocity  profiles  were  computed  such 
that  the  first  valid  bin  was  centered  at  7m  depth,  and 
velocities  where  estimated  for  each  4m  bin  when  75%  of  the 
acoustic  returns  were  good  within  each  1  minute  ensemble.  As 
the  raw  current  profiles  were  measured  in  the  ship's  reference 
frame  rotated  to  a  north/south  coordinate  system,  the  ship's 
velocity  must  be  removed  to  estimate  absolute  currents.  A 
commonly  used  procedure  is  to  choose  a  reference  layer  and 
subtract  its  velocity  from  all  the  other  bin  velocities.  This 
results  in  relative  velocity  profiles  free  of  the  ships 
velocity.  Similarly,  shear  profiles  will  not  be  contaminated 
by  the  ship's  velocity.  These  relative  velocity  data  are 
useful  for  studying  baroclinic  features,  like  the  shear  events 
mentioned  above;  however,  absolute  velocities  estimates  are 
needed  in  order  to  assess  barotropic  tides,  internal  wave 
fields,  mass  balances  over  the  full  closed  paths,  and  momentum 
and  vorticity  balances  at  specific  locations. 

The  major  problem  in  absolute  velocity  estimation  is 
to  obtain  an  accurate  estimate  of  the  ship's  velocity.  The 
ship's  mounted  GPS  navigation  system  has  a  50m  accuracy  for 
each  independent  position.  However,  by  analyzing  the  navigated 
ship's  velocity  spectral  characteristics  displayed  in  Figure 
II-B-2.1,  the  velocity  estimates  are  significantly  more 
accurate.  As  the  ship  was  usually  traveling  at  a  steady 
velocity  along  each  leg,  the  observed  high  frequency 
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variability  in  the  navigated  speed  was  dominated  by  the  GPS 
noise.  After  deglitching  the  data  using  a  threshold  value  such 
that  95%  of  the  data  is  assumed  good,  the  standard  deviation 
of  the  navigated  ship's  velocity  associated  with  high 
freguencies  (above  10  Hr'^)  is  approximately  O.lOm/s.  This  high 
frequency  velocity  variance  has  contributions  from  both  real 
physical  structure  and  the  ship's  positioning  errors.  Taking 
this  into  account,  the  GPS  ship  velocity  was  added  to  the  ADCP 
vertical  mean  velocities  to  form  an  absolute  velocity 
reference  layer,  after  the  reference  layer  was  filtered  by  a 
fourth  order,  zero-phase  Butterworth  filter  with  a  cut-off  at 
1/8  min',  to  filter  the  high  frequency  noise.  Another  source 
of  error,  which  must  be  taken  into  account  when  computing  the 
ADCP  (U,V)  components,  is  that  the  true  north  does  not 
correspond  exactly  to  the  direction  used  in  the  computations 
based  on  the  ship's  gyro  compass.  A  heading  correction  was 
computed  and  implemented  using  correlation  methods  (Cheriskin, 
1989) .  The  resultant  profiles  of  absolute  velocities  have  a 
95%  error  estimated  to  be  +/~0.04m/s,  based  on  the  RMS  value 
in  the  error  band  as  shown  in  Figure  II-B-2.1.  Due  to  the 
observed  highly  inhomogeneous  velocity  field,  some  of  the 
variance  in  the  error  band  is  probably  due  to  the  physical 
current  structure.  Therefore  this  error  value  slightly 
overestimates  the  true  uncertainty,  and  the  horizontal 
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velocity  gradients  are  slightly  smoothed  due  to  the  low  pass 
filter. 

With  data  in  this  gridded  form,  an  analysis  of  the 
complex  spatial  and  temporal  time  series  was  performed  to 
develop  the  hypotheses  which  will  be  discussed  in  the 
following  chapter. 
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III.  HYPOTHESIS  FORMULATION  BASED  ON  OBSERVATIONS 


A.  DATA  OVERVIEW 

This  section  presents  and  discusses  space-time  series  of 
the  velocity  and  temperature  fields  observed  during  the  MBAY 
experiment.  The  MBAY  data  are  difficult  to  represent  using 
traditional  methods,  so  maps  and  section  plots  corresponding 
to  snapshots,  or  time  averages,  of  14  different  data  fields 
have  been  used  to  present  the  highly  inhomogeneous  and  non¬ 
stationary  fields.  Several  complicated  plots  of  variables 
showing  examples  of  different  spatial  regions  and  times  are 
used  to  characterize  the  features,  which  dominated  the  upper 
ocean  during  the  measurement  period.  Due  to  the  nature  of  the 
display,  some  of  these  figures  show  features  that  are  not 
discussed  in  the  text  because  they  were  found  not  to  have  a 
significant  contribution  to  processes  relevant  to  the  present 
study . 

The  winds  were  moderate  with  magnitudes  (W)  below  5m/s 

throughout  the  entire  observation  period,  with  a  sharp  shift 

in  direction  at  hour  154.  The  resulting  wind  stress  was  less 
1  *2 

than  0.04m' s' Kg,  which  corresponds  to  an  unstratified  mixing 
depth  (h,)  of  the  order; 
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hg  =  lO'^xhf  *  10"^x^/£^x-p  »  20^7 


(8) 


where  Cq=1.4x10'^  is  the  drag  coefficient  and  f  is  the  Coriolis 
parameter . 

The  Ekman  solution  for  the  magnitude  of  the  surface 
velocity  (V^)  gives: 


Vo 


^  ySIlT 

Dgpf 


(9) 


where  r  is  the  wind  stress,  p  is  the  water  density  and  D,  is 
the  depth  of  frictional  influence  given  by: 


D 


(10) 


where  A,  is  the  eddy  viscosity.  Based  on  the  observed  winds 
and  for  an  homogeneous  ocean,  it  can  be  assumed  A^0.014m^s*^ 
or  smaller  (Pond  and  Pickard,  1983  -  p.l09) .  Applying  equation 
9  for  the  observed  conditions  gives  a  maximum  velocity  at  the 
ocean  surface  of  about  V^O.OSSms'^,  which  corresponds  to  a 
depth  of  frictional  influence  of  about  D,=5€m. 

However,  the  CTD  observations  showed  strong  stratification 
below  10m  (as  it  will  be  discussed  below) ,  suggesting  that  the 
wind  stress  was  distributed  over  a  shallower  depth.  Assuming 
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a  significant  mixed  layer  depth  of  h,=D.=10m,  the  magnitude  of 
the  surface  velocity  is  VQ=0.20ms'^  and  the  corresponding  slab 
current  magnitude  (i.e.,  velocity  integrated  over  the  depth  of 
frictional  influence)  is  about  0.04ms'\  Therefore,  taking  into 
account  the  7m  depth  first  bin  ADCP  current  observation  and 
the  0.04ms*’  accuracy  of  ADCP  current  observations,  wind  forced 
ocean  currents  are  not  likely  detected  in  the  MBAY  current 
data. 

Figure  III-A.l,  shows  the  time  series  of  the  winds  during 
the  experiment.  Vertical  lines  represent  the  starting  time  of 
each  path  around  the  observation  box  (total  of  27) ,  and  the 
hours  of  the  predicted  high  and  low  tides  are  marked  along  the 
x-axis.  These  figures  show  that  there  is  no  consistent 
spatial  variability  in  the  wind  between  consecutive  boxes, 
which  could  result  in  a  significant  wind  curl  capable  of 
forcing  vertical  velocities  at  the  mixed  layer  depth. 

Some  temporal  variability  can  be  observed,  particularly  in 
the  wind  direction,  at  hour  154.  Assuming  the  extreme  case  of 
wind  reversal  over  a  5km  length  scale,  the  maximum  resultant 
Ekman  pumping  would  force  vertical  velocities  at  the  mixed 
layer  depth  of  approximately. 


_ 0.03 _ 

5x10^x1024x0 . 87  xlO'* 


6xl0'®jms"^  (11) 
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producing  l.Sin  maximum  thermocline  displacements  if  sustained 
over  a  6  hour  period.  As  the  lateral  wind  gradients  could  not 
in  fact  be  detected,  W.  can  be  assumed  to  be  negligible  in 
this  study. 

The  objectively  analyzed,  vertically  averaged,  current 
data  temporally  averaged  over  the  full  3  days  of  the 
experiment  (Figure  III>A.2)  show  a  weak  anticyclonic  mean  flow 
of  approximately  0.05m.s'\  The  temporal  mean  net  velocity 
entering  the  closed  area  is  +0.0l8m/s  indicating  mean 
downwelling;  however,  this  value  is  comparable  to  the 
measurement  error.  Stronger  mean  currents  are  observed  along 
the  west  leg  going  northwestward,  while  weaker  mean  currents 
are  seen  along  the  east  leg,  but  in  the  opposite  direction.  In 
Figures  III-A.3  and  III-A.4,  the  spatial  structure  of  the 
temporal  mean  U  and  V  velocity  components  show  that  the 
stronger  currents  in  the  south-west  side  of  the  surveyed  area 
extend  from  25  to  190m,  reaching  a  maximum  of  about  0.25m/s  at 
70m  depth.  These  strong  mean  baroclinic  currents  are  not  seen 
on  the  north-east  side  of  the  surveyed  area. 

To  investigate  the  nature  and  possible  sources  of  these 
highly  inhomogeneous  mean  currents,  the  box  mean  divergence 
and  circulation  were  estimated.  These  estimates  are  very 
sensitive  to  the  ADCP  0.04m/s  RMS  velocity  error,  which  in 
this  case  is  dominated  by  the  positioning  system,  as  discussed 
in  the  previous  chapter.  Figure  III-A.5.a  shows  a  profile 
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time  series  of  the  depth  distribution  of  the  net  velocity 
entering  and  leaving  the  surveyed  area,  normalized  by  the 
Coriolis  parameter  and  divided  by  the  overall  area,  such  that 
the  units  correspond  to  vorticity  time  changes.  Magnitudes 
below  20xl0'^s'^  can  be  assumed  to  be  within  the  ADCP  noise 
level.  By  continuity,  the  vertically  averaged  time  series 
suggests  that  upward  and  downward  vertical  motions  are  present 
at  the  base  of  the  box,  approximately  at  tidal  periods, 
although  the  small  magnitude  of  the  vertical  mean  divergence 
is  not  above  the  noise  level  of  the  navigated  ADCP  profile 
data.  Nevertheless,  the  tidal  cycle  can  be  seen  in  the 
baroclinic  structure  of  the  estimates,  for  example  at  hour  150 
and  30m  depth,  negative  divergence  values  are  observed 
opposing  the  deeper  positive  values  along  the  same  profile 
with  magnitudes  well  above  the  noise  level.  At  hour  162 
positive  divergence  is  observed  between  20  and  130m  depth  with 
strong  negative  peaks  above  and  below,  and  at  hour  174  there 
is  a  strong  negative  divergence  near  the  surface  to  50m  depth 
and  positive  divergence  below.  This  divergence  maxima  should 
be  associated  with  strong  vertical  velocities  and  be  related 
to  time  increases  in  the  vorticity  magnitude.  In  Figure  III- 
A.5.b,  the  box  integrated  along-track  velocity  are  used  as  an 
indicator  of  the  average  vorticity  inside  the  surveyed  box, 
with  magnitudes  below  20xl0'^s'^  within  the  ADCP  noise  level. 
Strong  mean  vorticity  is  seen  in  the  profile  at  hour  150,  with 
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positive  maxima  at  20  and  170m  depths  and  a  negative  peak  at 
50m  depth.  The  vorticity  peaks  at  20  and  50m  depth  are 
consistent  with  the  negative  divergence  at  the  same  hour  at 
30m  depth.  On  the  other  hand,  vorticity  maxima  are  also 
observed  at  hour  134,  220m  depth  and  at  hour  172,  70m  depth, 
suggesting  that  the  layers  of  maximum  vorticity  are  shallowing 
in  time. 

If  the  flow  offshore  is  non-rotating  and  linear,  some 
local  source  must  be  forcing  the  observed  vorticity.  Based  on 
the  value  of  the  mixed  layer  slab  velocity  computed  above  and 
a  5Km  length  scale,  the  observed  light  and  steady  winds  until 
hour  154  can  produce  a  maximiim  divergence  by  Ekman  pumping  of 
6xl0'*s'^  (below  the  measurement  error) ,  so  the  observed  wind 
field  cannot  explain  the  changes  in  divergence  magnitude  at 
30m  depth,  and  the  observed  tidal  periodicity  of  the 
divergence  peaks.  Mechanisms  based  on  the  interaction  of 
internal  waves  with  tidal  flow  above  steep  topography  may 
explain  the  observed  vorticity  (and  divergence) ,  and  are 
analyzed  in  the  next  chapter. 

The  temporally  averaged  velocities  in  Figures  III-A.2/3 
have  a  strong  northward  component.  To  see  the  spatial  and 
temporal  structure  of  this  flow,  the  north-south  current 
component  is  displayed  in  Figure  III-A.6  at  four  equally 
spaced  points  along  the  west  leg.  These  figures  show  that  the 
strong  baroclinic  residual  velocities  are  associated  with 
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bursts  of  kinetic  energy  24  hours  apart,  most  clearly  seen  in 
the  profiles  at  location  '1560n'.  These  jet-like  events  have 
maxina  at  approximately  hour  142  (velocity  event  #1)  and  166 
(velocity  event  #2) ,  between  50  to  150m  and  30  to  100m  depth 
respectively.  They  are  unidirectional  (towards  northwest)  with 
0.30ms''  peak  velocities.  The  profiles  at  '2520m'  and  '3480m' 
show  the  first  event  with  two  local  maxima  in  depth,  one  at 
70m  and  the  other  at  150m;  however,  this  pattern  is  not  seen 
as  clearly  in  the  second  peak  at  hour  166.  The  profiles  at 
points  '600m'  and  '1560m'  suggest  that  similar  velocity  peaks 
occurred  a  few  hours  before  the  start  of  the  data  set  at  hour 
120  and  shortly  after  the  end  of  the  data  set,  although  this 
last  velocity  maximum  is  shallower  than  the  previous  features 
and  reaches  the  surface.  This  implies  that  the  strong  residual 
mean  velocities  in  the  west  pairt  of  the  box  are  a  direct 
consequence  of  baroclinic,  unidirectional,  jet-like  events 
occurring  every  24  hours,  which  shallow  with  time.  Such  jet 
events  did  not  occur  on  any  of  the  other  legs.  This  strong 
spatial  inhomogeneity  suggests  that  the  jets  could  be 
correlated  with  the  topography.  Unfortunately,  the  surveyed 
area  was  not  large  enough  to  determine  the  south-west  limit  of 
the  jets  and  to  positively  identify  the  maximum  currents. 
Inspection  of  the  four  panels  in  Figure  III-A.6  shows  that  the 
two  initial  events  at  hours  120  and  142,  and  the  last  event  at 
hour  186  are  stronger  at  locations  '600m'  and  '1560m',  however 
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the  event  at  hour  167  at  location  '600in'  is  strongly 
attenuated.  Coaparing  the  tiae  of  the  peaks  of  these  events 
and  the  predicted  tides  in  Monterey  Bay  (see  Figure  III-A.l), 
the  aaxiaua  velocities  occur  during  the  ebb,  corresponding  to 
the  aaxiaua  tidal  flow. 

Chapter  III-B  and  III-C  will  present  detailed  discussions 
of  the  observations  before  the  peak  of  the  event  at  hour  142 
(velocity  event  #1)  and  shortly  after  the  peak  of  the  event  at 
hour  166  (velocity  event  #2) . 

B.  VBLOCITY  BVBOT  «1 

This  section  discusses  the  following  features  of  the 
stratification  and  velocity  profile  tiae  series,  which 
characterize  the  physical  processes  in  the  surveyed  area  just 
prior  to  the  jet-like  feature  observed  at  hour  142: 

1  -  Theraocline  oscillations  with  tidal  periods  and 
aaplitudes  of  up  to  20a  occurring  at  the  north-west  comer  of 
the  surveyed  area. 

2  -  Strong  vertical  shear  occurring  below  the 
theraocline,  with  a  shallow  upper  layer  not  showing  the  saae 
fora  of  tidal  dependence  as  the  layer  below.  Unfortunately,  no 
temperature  (T)  and  salinity  (S)  data  were  available  at  the  peak 
of  this  event  or  deeper  than  80a. 

The  vertical  structure  of  temperature  around  the  seunple 
box  at  hours  124  (flood)  and  130  (ebb)  in  Figures  III-B.l  and 
III-B.2  show  theraocline  displacements  up  to  20a  during  the  6 
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hour  period  at  the  west  side  of  the  sampled  box  with 
wavelengths  larger  than  the  box  leg.  These  low  horizontal  mode 
displacements  have  larger  amplitudes  at  the  end  of  the  west 
leg  and  at  the  start  of  the  north  leg.  However,  their 
propagation  speed  or  wavelength  are  not  well  defined  due  to 
the  2  hour  sampling  time  around  the  box  perimeter.  Because 
there  are  no  obvious  mechanisms  to  cause  standing  waves  in  the 
area,  it  is  reasonable  to  assume  that  these  oscillations  are 
semi-diurnal,  propagating  internal  waves  with  wavelengths 
longer  than  lOKm  superimposed  on  the  mean  and  tidal  flows. 
Figure  III-B.3.b  displays  time  series  contours  of  the  V 
ccmponent  of  the  current  along  the  box  sides  at  19m  depth, 
within  the  upper  thermocline,  showing  a  positive  (northward) 
velocity  at  the  west  leg  between  hours  124  and  130. 
Furthermore,  along  the  west  leg,  velocity  changes  in  time  with 
a  semidiurnal  periodicity  are  consistent  with  the  period  of 
the  thermocline  oscillations  and  with  the  assumption  that 
these  oscillations  are  the  result  of  a  propagating  internal 
wave. 

This  type  of  thermocline  oscillation  is  not  as  clearly 
seen  in  the  other  parts  of  the  box;  however,  across  the  canyon 
(south  leg  at  hour  124,  south  and  north  leg  at  hour  130  in 
Figure  III-B.1/2) ,  the  spatial  distribution  of  other  groups  of 
smaller  scale  (-IKm)  internal  wave  isothenn  displacements  of 
about  10m  can  be  seen  superimposed  on  the  larger  scale 
features.  The  measurement  time-space  resolution  and  time 
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series  lengths  are  not  adequate  to  estimate  their  phase  or 
group  speeds. 

The  temperature  profiles  displayed  in  Figures  III-B.l  and 
III-B.2  show  the  upper  ocean  is  stratified  to  the  surface  so 
that  no  surface  mixed  layer  is  seen  over  much  of  the  area.  In 
fact,  strong  upper  layer  spatial  inhomogeneities  in  the 
temperature  profiles  occur  over  the  south  and  west  legs.  This 
stratified  surface  layer  is  consistent  with  upwelling 
conditions  and  the  observed  light  winds  shown  in  Figure  III- 

A. l,  which  persisted  throughout  the  measurements. 

The  current  time  series  contour  maps  of  Figures  III- 

B. 4.a/b  and  III-B.5.a/b  summarize  the  spatial  and  temporal 
currents  at  11m  and  51m  depths.  Strong  differences  are  evident 
between  them.  Before  hour  156,  the  U  component  of  the  current 
displayed  in  Figure  III-B.4.a  suggests  a  mean  value  of  about 
5cm/s  at  11m  depth  superimposed  on  the  tidal  cycle.  The  tidal 
current  is  most  clearly  seen  in  the  V  component  plot  in  Figure 
III-B.4.b.  For  the  same  period.  Figure  III-B.5.b  shows  strong 
diurnal  northward  jets  at  51m  depth  along  the  west  leg 
superimposed  with  regular  tidal  cycles,  seen  in  the  U 
component  plots  of  Figure  IIl-B.5.a.  The  different  behavior  of 
the  upper  and  lower  layers  before  hour  156  suggests  that 
different  mechanisms  are  driving  the  residual  currents,  though 
both  appear  to  be  related  to  the  tides. 


C.  VBLOCITY  BVBMT  *2 

This  section  discusses  the  following  features: 

1  -  A  velocity  peak  at  hour  167  with  different 
characteristics  from  those  discussed  in  the  previous  section. 

2  >  An  internal  mixed  layer  between  20m  and  50m  depth. 

3  -  The  development  of  this  internal  mixed  layer  at 
hour  160,  prior  to  the  velocity  peak  of  the  event  and 
following  the  wind  shift  at  hour  154.  There  were  no 
temperature  and  salinity  data  available  between  hours  134  and 
166,  or  below  80m  depth. 

The  large  scale  internal  wave  displacements  of  the 

thermocline  described  in  the  previous  section  are  not  seen  in 

the  vett:ical  section  snapshots  of  temperature  in  Figures  III- 

C.l  at  hour  168  and  III-C.2  at  hour  174.  However,  an  internal 

homogenized  layer  is  evident  in  the  west  and  north  legs  at 

hour  168,  between  15  and  50m  depth  (Figure  III-C.l).  This 

internal  mixed  layer  is  well  developed  at  hour  162,  prior  to 

the  peak  of  the  event  at  hour  166  (see  Figure  III-C.3) .  Later, 

by  hour  174,  this  layer  has  restratified  as  the  velocities 

decreased  (see  Figure  1II-C.2  west  and  north  legs). 

» 

Below  and  above  this  mixed  layer,  strong  vertical 
temperature  gradients  are  seen  in  the  isotherm  structure  of 
Figures  III-C.l  south  and  west  legs,  and  Figure  III-C.2  west 
and  north  legs.  These  figures  also  suggest  that  smaller  scale 
internal  waves  might  be  present,  but  their  signatures  are  not 
as  obvious  as  those  shown  in  the  previous  section. 


A  patch  of  warm  water  shown  in  Figure  III-C.!,  (west  leg 
at  position  lOOOm)  ,  displaced  the  12*’c  isotherm  downward  by 
almost  40m.  During  this  event  at  hour  168,  the  velocity  jet 
following  a  24  hour  cycle  along  the  west  leg  did  not  achieve 
the  same  magnitude  as  the  previous  jets  (see  Figure  III-A.6- 
600m/2520  and  3480  panels) . 

The  high  magnitude  8  meter-scale  vertical  shear  between 
30m  and  the  surface,  displayed  in  Figure  III-C.3  at  hour  156 
in  the  profiles  at  point  '600m'  in  the  west  leg,  is  likely  to 
have  triggered  the  internal  mixing  event.  This  mixed  patch  was 
then  seen  as  an  internal  mixed  layer  in  the  west  leg  when  the 
T/S  Tow-Yo  measurements  were  resumed  at  hour  168,  as  displayed 
in  Figure  III-C.l. 

Figure  III-B.3.b  displays  time  series  contours  of  the  V 
component  of  the  current  along  the  box  sides  at  19  meters 
depth,  within  the  upper  thermocline  and  near  the  depth  of  the 
mixing  event  observed  in  Figures  III-C.1/2.  The  contours  show 
that  at  approximately  hour  156  there  was  a  reversal  (from 
northward  to  southward  flow)  of  the  north-south  velocity 
component.  This  temporal  change  is  in  agreement  with  the  upper 
layer  velocity  at  11m  depth,  as  displayed  in  Figure  III-B.4.b 
at  hour  156,  but  is  not  observed  in  the  deeper  layers  at  51m 
depth  (Figure  III-B.5.b)  suggesting  that  the  mixing  event  was 
triggered  at  this  time  due  to  a  shear  instability  process 
related  to  an  enhanced  vertical  shear.  These  observations 
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suggest  that  the  observed  internal  mixed  layer  at  hour  168  was 
associated  with  an  instability  mechanism  at  the  thermocline 
related  to  the  presence  of  the  strong  jet-like  events. 
However,  the  data  do  not  allow  any  conclusions  about  the  exact 
triggering  mechanism. 

D.  HYVOTSESX8  TORMULXTXOM 

Based  on  the  data  discussed  in  the  previous  section,  the 
measured  weak  winds  cannot  force  the  observed  highly 
Inhomogeneous  and  strong  current  field,  high  vertical  current 
shear,  or  the  observed  semidiurnal  thermocline  oscillations. 

On  other  hand,  the  fact  that  the  measurement  site  was 
close  to  the  mouth  of  the  Monterey  canyon  suggest  that  canyon 
dynamics  might  be  relevant.  The  effects  of  geostrophic  flow 
over  a  canyon  has  been  studied  by  Klinck  (1988,  1989).  He 
showed  that  narrow  canyons,  defined  as  those  whose  width  is 
small  compared  to  the  internal  Rossby  radius  of  deformation 
(-15Km  for  Monterey  canyon)  (Koehler,  1990) ,  will  act  as  a 
barrier  to  barotropic  transverse  flow.  Using  the  results 
presented  by  Klinck  (1988)  and  the  data  of  Koehler  (1990) ,  the 
time  for  geostrophic  adjustment  of  the  flow  crossing  the 
Monterey  canyon  can  be  estimated  to  be  on  the  order  of  1.5  hr, 
exciting  waves  over  the  canyon  with  frequencies  4  and  9  times 
the  inertial  frequency,  and  radiating  waves  from  each  side 
with  frequencies  5  to  13  times  the  inertial  frequency.  The 
effects  of  stratification  were  discussed  by  Klinck  (1989) .  For 
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narrow  canyons  like  Monterey  canyon,  he  predicts  slower 
velocities  over  the  canyon  with  a  cyclonic  circulation  within 
the  canyon  and  an  anticyclonic  circulation  above.  In  contrast 
to  the  case  of  along-canyon  flow  (like  the  tidal  flow) ,  Klinck 
(1988,1989)  suggests  that  the  canyon  will  always  affect  the 
external  pressure  gradients,  generating  residual  flows  over 
the  shelf.  However,  for  the  area  covered  in  our  survey,  the 
cross'canyon  geostrophic  flow  is  expected  to  be  much  less 
energetic  than  the  tidal  flow,  based  on  previous  work  in  the 
area  summarized  in  Breaker  et  al.,  1989.  This  is  confirmed  by 
the  unpublished  current  meter  buoy  data,  kindly  made  available 
by  C.  H.  Pllskaln  (Pillsbury,  et  al.,  1992)  (see  Figures  III- 
D-l.a-d),  which  show  strong  tidal  currents  aligned  along  the 
canyon  axis.  As  a  result,  the  effects  of  geostrophic  flow 
across  the  canyon  on  the  generation  of  waves  and  residual 
flows  is  not  likely  to  explain  the  MBAY  observations. 

Therefore,  based  on  the  discussion  of  the  data  presented 
in  the  previous  section,  the  following  hypotheses  are 
established: 

-  Internal  tides  are  being  observed  at  the  box  survey 
site. 

-  The  unidirectional,  baroclinic,  jet-like  features 
are  the  result  of  rectifications  of  internal  tides 
propagating  below  the  thermocline 
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-  The  internal  well-mixed  layer  discussed  in  the 
previous  section  is  a  result  of  an  instability 


mechanism  of  the  internal  tide  at  the  thermocline. 

A  detailed  discussion  of  these  hypotheses  follows. 

1.  Tidally  Generated  Internal  Waves 

The  barotropic  tide  interacting  with  the  shaJllowing 
continental  shelf  topography  forces  vertical  velocities  as  the 
horizontal  tidal  current  magnitude  increases.  The  pressure 
gradients  established  by  the  response  of  the  stratified  water 
column  to  these  vertical  velocities  will  force  internal  waves 
with  tidal  periods  (internal  tides).  These  internal  waves 
propagate  away  from  the  shelf  break,  where  the  topographically 
generated  pressure  gradient  forcing  is  maximum  (Baines,  1982; 
Huthnance,  1989).  The  observed  long  wavelength  thermocline 
oscillations  discussed  in  IIl-B  are  consistent  with  a 
baroclinic  interfacial  mode  within  a  2-layer  like  upper  ocean 
structure . 

A.  Review 

Observations  of  internal  wave  perturbations 
propagating  along  the  axis  of  the  Monterey  canyon  have  been 
reported  by  Shepard  et  al.  (1974).  They  report  phase 
velocities  of  about  0.25ms'^  for  deeper  water  and  of  about 
0.38ms'’  in  shallower  water.  Internal  tides  over  the  California 
continental  shelf  have  been  observed  at  several  locations 
(Huthnance,  1989).  Shea  and  Broenkow  (1982)  observed  large 
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tidal ly  related  isotherm  displacements  at  the  head  of  the 
Monterey  canyon.  Their  observations  suggest  that  strong 
internal  tides  should  exist  in  the  area,  but  their  data  were 
not  sufficient  to  allow  a  full  regional  characterization. 

Baines  (1983)  used  a  laboratory  experiment  to  study 
tidal  motions  over  submarine  canyons.  He  described  a  process 
by  which  the  internal  tide  amplitudes  can  be  enhanced  inside 
narrow  canyons  (defined  as  canyons  whose  width  is  much  smaller 
than  both  the  tidal  wavelength  (~200Km)  and  the  barotropic 
Rossby  radius  of  deformation) ,  so  that  the  canyon  will  have 
negligible  effect  on  the  pressure  gradients  established  by  the 
barotropic  motion.  Therefore,  the  tidal  motions  inside 
topographical  features  like  the  Monterey  canyon  will  be  forced 
by  the  same  pressure  gradients  that  exist  over  the  continental 
shelf,  and  by  a  common  pressure  gradient  established  at  the 
canyon  head,  leading  to  an  enhanced  transport  inside  the 
canyon  (Huthnance,  1989).  This  type  of  pressure  gradient 
supports  the  conceptual  model  proposed  by  Shea  and  Broenkow 
(1982)  of  cold  deep  canyon  water  being  lifted  by  the  tides  and 
intruding  onto  the  shelf  by  the  established  horizontal 
pressure  gradients.  At  the  MBAY  measurement  site  the 
topographic  slope  is  much  larger  than  the  internal  wave  slope, 
as  will  be  shown  in  the  next  chapter.  For  this  situation, 
Baines  (1983)  predicted  an  enhanced  internal  tide  up-canyon 
and  near  the  bottom,  verified  by  observations  of  strong 
internal  tides  by  Shea  and  Broenkow  (1982)  at  the  canyon  head. 
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He  also  predicts  that  the  tidal  flow  should  be  strongly 
baroclinic.  For  this  case,  Huthnance  (1989)  also  suggests  that 
the  tidal  flow  will  show  a  complex  spatial  dependence.  Based 
on  the  theory  by  Grimshaw  et  al.  (1985),  reflections  at  the 
open  mouth  can  contribute  to  this  intensification,  especially 
for  steep  slopes  and  narrow  canyons. 

Because  the  MBAY  survey  area  is  located  close  to 
the  canyon  mouth  (but  not  within  the  canyon)  ,  it  can  be 
assumed  that  the  internal  tide  dynamics  will  be  very  similar 
to  those  over  the  open  shelf,  and  the  effects  of  the  canyon 
can  be  viewed  as  having  higher  order  effects.  Current  meter 
buoy  data  (Pillsbury,  et  al.,  1992)  displayed  in  Figures  III- 
D-l.a-d,  show  the  tidal  flow  at  500m  depth  close  to  the 
surveyed  area  (as  shown  in  the  map  in  Figure  II-A.l)  to  be 
largely  aligned  along  the  canyon  axis.  Significant  isotherm 
displacements  with  tidal  periods  are  clearly  seen  in  the 
temperature  data,  confirming  the  importance  of  tidal  dynamics 
in  the  area. 

B.  Problem 

To  identify  the  sources  for  internal  waves  and  to 
characterize  the  mesoscale  and  tidal  circulation  in  the  area, 
simultaneous  observations  would  be  needed  at  several  locations 
along  the  canyon,  along  the  shelf,  and  on  the  continental 
slope.  Nevertheless,  the  significant  along  canyon  tidal  flow 
identified  in  both  these  data  and  in  the  current  meter  mooring 
data  allows  some  simplifying  asstunptions  to  be  made  in 
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isolating  the  tides  as  a  primarly  source  of  internal  waves  in 
the  area,  other  mechanisms  will  be  assumed  to  modulate  the 
tidal  dynamics  on  a  longer  term  (at  least  several  days)  or  as 
higher  order  much  less  energetic  effects. 

Other  potential  problems  in  the  analysis  of 
internal  tidal  forcing  is  the  possibility  of  multiple 
reflection  effects  inside  the  canyon  and  modification  of  the 
internal  wave  field  by  pressure  gradients  established  by 
residual  flows.  However,  because  our  measurement  box  is  at  the 
canyon  mouth,  the  assumption  that  the  observed  dynamics  are 
similar  to  those  on  the  continental  slope  is  justified  to  a 
first  approximation  and  the  effects  of  local  pressure  fields 
can  be  regarded  as  an  higher  order  correction,  superimposed  on 
the  larger  scale  response. 

The  characterization  of  the  barotropic  tides  in  the 
area  will  also  require  some  special  care.  Unfortunately,  no 
complete  regional  study  of  the  tides  exist  for  the  area, 
c.  Approach 

To  evaluate  the  internal  tide  forcing  in  the  area, 
the  linear  model  developed  by  Baines  (1982)  is  used  in  the 
next  chapter  to  estimate  the  magnitude  of  the  forced  internal 
wave  field  on  a  regular,  straight,  continental  slope  with 
representative  model  parameters  for  the  survey  area. 

To  better  simulate  the  observed  data  and  to  take 
into  account  some  of  the  enhancements  which  might  occur  due  to 
the  presence  of  the  canyon,  the  barotropic  tidal  flux  will  be 
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assumed  to  have  0.05m/s  amplitudes  at  1000m  depth,  increasing 
as  depth  decreases,  and  generating  vertical  velocities  at  the 
bottom  boundary  by  continuity.  The  resulting  vertical 
velocities  decrease  linearly  with  height  towards  zero 
magi  .de  at  the  sea  surface. 

2 .  Tidal  Reotif loatlon 

Non-linear  effects  of  the  barotropic  tides  interacting 
with  topography  and  the  internal  tides  result  in  spatially 
inhomogeneous  residual  currents.  Therefore,  it  is  hypothesized 
that  the  unidirectional,  baroclinic,  jet  like  features 
observed  on  the  west  leg  of  the  surveyed  box  are  caused  by 
rectification  of  the  internal  tide  due  to  non-linear  effects. 
It  is  further  hypothesized  that  the  nonlinearities  are 
dominated  by  internal  tides  interacting  with  the 
topographically  enhanced  barotropic  tidal  flow  (assumed  to  be 
a  free  shallow-water  gravity  wave)  and  the  water  column 
stratification.  These  hypotheses  will  be  justified  in  the 
following  sections,  and  tested  in  a  model  described  in  Chapter 
IV. 

A.  Review 

Tides  are  a  periodic  phenomena,  which  ideally  have 
zero  mean  velocities  over  long  periods.  However,  due  either  to 
ocean  floor  topography,  frictional  effects,  or  to  the  coastal 
shape,  there  are  processes  that  can  rectify  tidal  currents 
resulting  in  spatially  inhomogeneous  mean  currents.  Tidal 
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r«c±ification  has  been  shown  to  be  due  to  non-linear  effects 
and  friction  in  a  homogeneous  ocean  (Robinson,  1985;  Mass  et 
al.,  1987;  Garreau  et  al.,  1992),  and  to  non-linear  effects  on 
a  baroclinic  tidal  flow  associated  with  the  internal  tides 
(Maze,  1987;  Willmott,  1987;  Pichon  et  al.,  1990;  Mass  et  al., 
1989a, b) . 

Following  Robinson  (1985),  rectified  flows  are 
typically  one  or  two  orders  of  magnitude  less  than  the  tidal 
current  velocities,  requiring  special  care  in  observation  and 
modeling.  Robinson  Identifies  the  non-linear  interactions  as 
the  main  process  by  which  residual  vorticity  can  occur,  and 
presents  four  vorticity  generation  mechanisms  which  play  an 
impoxrbant  role  in  coastal  dynamics,  as  shown  in  the  Figure 
III-D.2.a; 

-  potential  vorticity  conservation  (column 
stretching  and  squeezing) ; 

-  differential  friction  due  to  horizontal  shear; 

-  differential  friction  due  to  lateral  depth 
gradients ; 

-  asymmetrical  depth  profile  changes. 

For  surface  displacements  small  compared  to  depth, 
small  residual  flows  result  (compared  to  the  tidal  currents) 
with  a  relative  vorticity  much  larger  than  the  planetary 
vorticity  (assumptions  likely  to  be  valid  for  the  MBAY  data 
considering  the  vorticity  estimates  of  Figure  III-A.6). 
Zimmermann  (1978,  1981)  used  a  simplified  potential  voirticity 
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equation  to  conclude  that  residuals  are  determined  by  a 
balance  between  generation  by  tidal  advection  and  dissipation 
by  turbulent  eddy  diffusion.  He  found  that  maximum  residual 
vorticity  occurs  over  topographic  features  with  length  scales 
comparable  to  the  tidal  excursion.  For  the  case  of  the 
continental  slope  near  Monterey,  the  typical  topographical 
length  scale  of  about  5Km  is  comparable  to  a  tidal  excursion 
length  of  about  3.5Km  for  lOcm/s  tidal  currents,  suggesting 
that  these  mechanisms  might  be  relevant  in  the  area.  The  tidal 
excursion  length  (Lf)  is  defined  as 


Ly  =  2 


^tiae(x) 

<d 


(12) 


where  is  the  tidal  frequency.  Using  similar  assumptions, 
Robinson  (1981)  observed  using  simplified  niimerical  solutions 
in  a  quasi-lagrangian  sense  the  following  features: 

-  vorticity  effects  cannot  penetrate  further  than 
the  tidal  excursion; 

>  vorticity  can  propagate  normal  to  the  tidal 
streamlines  only  by  eddy  diffusion  processes; 

-  because  vorticity  is  dissipated  within  one 
tidal  time  scale,  it  cannot  be  advected  further 
than  the  length  of  the  tidal  excursion,  so  the 
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areas  of  influence  of  topographic  features  is 
limited; 

-  tidal  advection  of  a  gradient  of  tidal 
vorticity  is  the  mechanism  by  which  residual 
vorticity  is  generated. 

Also  by  running  simplified  numerical  solutions  in  a  Eulerian 
sense,  including  only  the  stretching  mechanisms  terms,  he 
found  that: 

-  positive  mean  vorticity  is  generated  in  deeper 
water  and  negative  mean  vorticity  is  generated  at 
shallower  depths; 

-  the  maximiun  generation  of  residual  vorticity 
occurs  when  the  topographic  length  scale  is 
similar  to  the  local  tidal  excursion. 

He  estimated  the  maximum  residual  vorticity  ((q)  ^  shallow 
sea  due  to  column  stretching  to  be 


0.45 


£V  dh 
ah  ds 


(13) 


In  the  same  way,  considering  only  bottom  friction,  he  found 
that  the  maximum  residual  vorticity  in  shallow  seas  to  be 


0.4 


c^v  a 

6)  dn 


(14) 
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^•rm  'V  is  the  tidal  stream  amplitude  in  the  's' 
(streamline)  direction,  'n'  is  the  direction  normal  to  's', 
**Cb"  is  the  bottom  friction  coefficient,  'f  is  the  planetary 
vorticity,  and  'h*  is  the  water  depth.  Based  on  these  results, 
residual  vorticity  patterns  can  be  dra%m  for  a  specific 
location  and  the  mean  residual  flow  estimated. 

In  the  case  of  Monterey  Bay,  depths  range  from  1000m 
to  100m  over  distances  of  5000m,  corresponding  to  very  large 
slopes  of  approximately  15%  and  100%  depth  changes  within  one 
tidal  excursion  (~3.5Km).  However,  despite  this  strong 
topography,  barotropic  tidal  currents  at  these  depths  ranges 
are  expected  to  be  smaller  than  0.50ms'\  so  that  column 
stretching  will  be  the  mechanism  likely  to  generate  tidal 
residual  flow.  The  rectified  flow  will  be  maximized  for 
topographic  lengths  ec[ual  to  the  tidal  excursion  i.e.,  As«V/». 
For  this  case  and  for  typical  values  in  the  continental  slope 
close  to  Monterey,  the  maximum  residual  vorticity  will  be 
bounded  by  (,«0.45f  (Ah/h)  <  3.€xlo'*s'\  corresponding  to  maximum 
residual  flows  smaller  than  v«((0.x«f)/4  ••  O.OSms*^  for  shallow 
waters  and  ▼••O.OOSms'^  for  the  deep  water.  Figure  III-A.5  shows 
a  residual  mean  vorticity  of  about  3xl0*^s'^  in  the  surveyed 
area,  indicating  a  good  agreement  with  this  theory.  However, 
the  observed  residual  vorticity  has  a  strong  depth  dependence 
indicating  that  the  barotropic  assumption  is  not  fully  valid. 
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This  slow  rectified  flow  will  result  from  northward 
currents  along  the  continental  slope  with  shallow  water  on  the 
right.  This  simple  analysis  suggests  that  barotropic  tidal 
rectification  may  play  some  role  in  the  dynamics  of  Monterey 
Bay.  However,  the  small  magnitudes  of  the  rectified  flows  fall 
within  0.02-0.04ms'^  noise  floor  of  our  observations. 

More  accurate  theoretical  and  numerical  models  for 
barotropic  tidal  rectification  have  been  developed  by  Mass  et 
al.  (1987)  and  Garreau  (1992)  among  others.  They  better 
simulate  the  effects  of  stronger  non-linearities  that  are 
expected  to  occur  over  steep  slopes.  Garreau  (1992)  defines  a 
non-dimensional  parameter  indicating  non-linearity  as  Lr/Iro 
where 

L,^  =  il(x)/-g  (15) 

A  typical  value  of  L|-o  for  the  continental  shelf  close  to 
Monterey  Bay  is  35Km,  so  Lr/Lfo  order  of  0.1, 
comparable  to  values  in  the  Celtic  Sea  where  strong  and 
spatially  inhomogeneous  mean  currents  of  the  order  of  O.OSms*' 
have  been  observed.  The  Celtic  sea  residual  flows  have  been 
attributed  to  tidal  rectification  processes  and  this  site  also 
has  internal  tides  predicted  (Baines,  1982)  and  observed 
(Huthnance,  1989) ,  although  different  length  scales  apply. 
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Non-linear  layer  nodels  have  been  developed  by  Maze 
(1987)  and  Willmott  (1987)  to  study  the  effects  of  advection 
by  the  barotropic  tides.  Significant  residual  currents  are 
predicted  in  the  upper  layer  changing  within  the  length  scales 
of  topography  (L^)  and  the  internal  wave  scale  (depth/c) , 
where  'c*  corresponds  to  the  slope  of  internal  waves  for  a 
typical  stratification.  It  will  be  shown  in  the  next  chapter 
that  this  last  scale  is  0(3Km)  for  the  data  set  under 
discussion.  The  simplest  analysis  for  weak  non-linearities 
predicts  typical  mean  currents  to  be  on  the  order  of 
(Huthnance,  1989) : 


where  'A*  corresponds  to  the  internal  wave  current  amplitudes 
and  'L*  the  significant  length  scale.  It  will  be  shown  in  the 
next  chapter  that  current  oscillations  are  predicted  with 
amplitudes  of  about  O.lms'^  for  the  interfacial  mode  in  the 
MBAY  data  set.  This  amplitude  will  result  in  residual 
velocities  of  0.02ms*^  in  the  upper  layer  and  0.002ms'^  in  the 
lower  layer  of  an  idealized  2  layer  system.  The  velocities  of 
the  forced  internal  waves  in  the  stratified  continuum  below 
the  thermocline  will  be  estimated  in  the  next  chapter  to  be 
approximately  0.15ms'\  resulting  in  residuals  of  about 
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0.09a8'\  Therefore,  advection  by  the  tidal  flow  appears  as  a 
significant  source  for  baroclinic  tidal  rectification. 

The  non-linear  modal  approach  for  a  stratified 
ocean  by  Mass  and  Zimmemann  (1989a, b)  showed  that,  though 
there  were  many  features  in  common  with  the  homogeneous  case 
developed  above,  a  bottom  residual  current  intensification  can 
be  also  observed.  They  found  a  depth  dependent  relative 
vorticity  produced  by  vortex  stretching  of  free  internal  modes 
causing  bottom  enhanced  residual  vorticity.  In  fact,  the 
vertical  structure  of  vorticity  was  found  by  them  to  be 
exponential  rather  than  sinusoidal.  They  also  found  evidence 
for  bottom  intensified  cross-isobath  circulation  due  to 
baroclinicity  and  along-slope  baroclinic  currents  extending 
for  greater  distances  than  for  the  barotropic  counterpart. 
Their  results  are  supported  by  observations  from  the 
continental  slope  near  northwest  Africa  discussed  by  Gordon 
(1980)  and  Huthnance  et  al.  (1982)  where  the  internal  tide  was 
observed  to  be  bottom  intensified.  Although  this  mechanism 
supports  bottom  current  intensification,  it  does  not  predict 
the  current  profiles  observed  in  the  MBAY  experiment  where  the 
stronger  currents  were  observed  between  50  and  150m  depths. 
Also,  the  buoy  data  at  500m,  described  above,  does  not  show 
any  evidence  of  intensification  with  depth  as  would  be 
required  if  Hass  and  Zimmermann  mechanisms  were  to  be  relevant 
to  the  MBAY  observations.  Thus,  although  these  measurements  do 
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not  allow  any  conclusion  about  the  existance  of  bottom 
Intensified  currents,  this  mechanism  cannot  explain  the  MBAY 
observations . 

B.  Problem 

None  of  the  theories  mentioned  above  seem  to 
account  for  the  occurrence  of  the  observed,  strong, 
unidirectional,  diurnal,  baroclinic  jets,  although  the  non¬ 
linear  interaction  of  finite  amplitude  internal  waves 
generated  at  the  thermocline  with  the  barotropic  tidal  flow 
appears  to  be  a  likely  mechanism.  Because  the  currents  of  the 
observed  jets  are  approximately  0.30ms'',  some  amplification 
mechanism  must  also  exist.  As  these  jets  did  not  occur  as  a 
bottom  current  enhancement  (they  decay  below  150m  depth) 
mechanisms  other  than  the  baroclinic  vortex  stretching 
discussed  by  Mass  et  al.  (1989)  are  present. 

C.  Approach 

Kunze  (1985)  showed  the  importance  of  including 
straining  together  with  advection  by  the  mean  flow  when 
analyzing  the  propagation  of  near-inertial  internal  waves.  A 
ray  tracing  model  based  on  this  approach  is  developed  in  the 
next  chapter  to  study  the  evolution  of  baroclinic  internal 
waves  generated  close  to  the  thermocline  by  the  barotropic 
tides,  using  an  idealized  continental  slope  with  values 
representative  of  those  either  side  of  Monterey  canyon.  The 
extension  to  Kunze 's  work  is  the  inclusion  of  vertical 
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boundaries  and  the  time  dependent  component  in  the  mean  flow, 
which  in  our  case  is  assumed  to  be  the  barotroplc  tide.  Also 
a  second  order  result  is  calculated  by  computing  the  wave- 
action  along  the  rays,  defined  as  the  ratio  of  the  wave  energy 
to  the  Euler ian  frequency.  Wave-action  was  proved  to  be  a 
Riemann  invariant  (see  Whitham,  1974)  for  internal  waves  by 
Bretherthon  et  al.  (1969),  and  can  be  used  to  estimate  the 
wave  amplitude  along  characteristics. 

3.  Znteimal  Tide  Instabilities 

The  last  hypothesis  is  that  the  vertical  shear 
supplied  by  baroclinic  rectifications,  both  on  the  interfacial 
thermocline  mode  and  on  the  internal  waves  that  are 
propagating  in  the  stratified  continuum  below,  will  supply 
sufficient  energy  for  shear  instabilities  to  occur  within  the 
thermocline.  Several  examples  of  these  strong  internal  mixing 
events  were  observed  in  the  MBAY  measurements  (see  Figure  III- 
C.l) . 

A.  Review 

Internal  waves  are  oscillations  about  an 
equilibrium  state  within  the  interior  of  a  stratified  fluid. 
Their  main  restoring  force  is  buoyancy,  and  depending  on  their 
length  scales,  rotational  effects  in  an  f-plane  will  be 
significant  for  near-inertial  frequency  internal  waves  like 
internal  tides.  The  presence  of  lateral  boundaries  and 
topography  may  introduce  extra  constraints  (LeBlciid  and  Mysak, 
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1972;  Gill,  1982).  Internal  waves  in  the  open  ocean  have  been 
extensively  observed  and  studied  in  the  past  twenty  years  (see 
for  example  Levine,  1983;  Muller,  et  al.,  1991). 

Based  on  the  current  knowledge  about  internal 
waves.  Garret  and  Munk  (1975)  developed  an  analytical,  random 
phase,  linear  model  (GM75)  describing  a  background  continuum 
of  wave  energy  that  has  been  successfully  compared  with 
oceanic  observations,  and  is  now  a  widely  accepted  standard 
for  comparison  and  analysis.  Deviations  from  this  background 
continuum  can  usually  be  associated  with  particular  generation 
mechanisms.  The  evolution  towards  a  GM75  spectral  shape  is 
explained  by  three  classes  of  non-linear  resonant  interacting 
triads:  elastic  scattering  (forced  vertical  symmetry), 
parametric  sub-harmonic  instability  (transfer  of  energy 
towards  smaller  scales  at  near-inertial  frec[uencies)  and 
induced  diffusion  (diffuses  wave  action  in  wavenumber  space 
forcing  the  high  wavenumber  spectrum  to  an  equilibrium  shape) 
(McComas,  et  al.,  1977;  Phillips,  1960). 

The  frequency  band  for  internal  waves  has  a  range 
limited  between  the  inertial  and  the  Brunt-Vaisaala 
frequencies,  although  near- inertial  waves  with  frequencies 
lower  than  inertial  (Kunze,  1985),  and  waves  with  frequencies 
above  N  (Lamb,  et  al.,  1992)  are  possible  due  to  wave-mean 
flow  interactions  and  horizontal  inhomogeneities. 

Non-linear  interactions  of  internal  waves  have  been 
identified  as  possible  sources  of  mixing  within  the  ocean 
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interior  (Gregg,  et  al.,  1986;  Caldwell,  1987;  New,  1988; 
Kunze,  et  al.,  1992). 


In  coastal  waters  strong  departures  from  the  GM75 
spectrum  might  be  expected  due  to  the  generation  of  internal 
tides,  the  effects  of  topography,  the  presence  of  a  lateral 
boundary  and  complicated  wind  regimes,  all  of  which  are 
thought  to  change  the  wave  field  and  induce  strongly  non¬ 
linear  conditions. 

The  following  analysis  follows  Leblond  and  Mysak 
(1978,  sections  42  and  43)  in  the  discussion  of  the  basic 
theory  for  stability  of  internal  waves,  neglecting  rotation. 
The  Taylor-Goldstein  equation,  derived  from  the  two- 
dimensional,  linearized  perturbation  equations  (p.  398), 
describes  internal  waves  superimposed  on  a  mean  vertical 
shear: 

[Po(U  -  Ov/]'  -  (PqU'w)'  -  [ -  +  Pok^U  -  C)3w  =  0 

u  -  c  ” 


where  ’Po<s)  ’  is  the  non-perturbed  density,  '0(z>  '  is  the  mean 
current,  'C  is  the  phase  velocity,  'w*  is  the  vertical 
velocity,  'g'  is  the  gravity  acceleration,  'k'  is  the 
horizontal  wavenumber,  and  primes  represent  vertical 
derivatives.  This  ec[uation  can  be  used  to  describe  the  onset 
of  shear  flow  Instabilities  as  predicted  by  a  linear  analysis. 
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For  a  two  layer  ocean,  with  U=UfS/h  (h  being  the  layer  depth) 
in  the  upper  layer  and  no  mean  flow  in  the  deeper  layer, 
equation  (17)  reduces  to: 

w"  -  ]cH/  =  0  (18) 


in  each  layer,  providing  U-C  *  0.  This  yields  solutions  of  the 
type: 


Ae'“  +Be'*'  ,  0  <  2  <  H 
De*'  ,  -00  <  z  <  0 


(18) 


where  A,  B  and  D  are  constants  determined  by  the  matching 
conditions  at  the  interface: 


[ 


=  0 


[Pot(U  '  C)w'  -  U'w  -  gw/(U  -  C)]]  =  0 


(20) 


The  kinematic  and  pressure  conditions  at  the  surface  combine 
to  give: 

(U  -  C)  V  -  (U  -  C)u'w  -  gw  =  0  (21) 

For  a  non  trivial  solution,  they  arrived  at  the  following 
eigenvalue  equation  for  ' c ' : 


9 


[  (U,  -  C)^k  -  g)  (A(c^k  -  g)  *  (U,  -  O^k  +  g)  + 

+  e-^[(U,  -  o’^k  +  g3[A(c^k  -  g)  -  [  (U,  -  C)^k  -  g]]  =  0 


(22) 


where  A  >  1  is  the  ratio  between  the  layer  densities.  This 
equation  has  four  real  roots,  or  two  real  roots  and  a  complex 
conjugate  pair  of  roots.  In  the  latter  case,  one  of  the 
complex  roots  generates  a  time  growing  (unstable)  solution, 
such  that  the  flow  becomes  unstable.  Therefore,  the 
st2d>ilizing  effect  of  stratification  is  not  always  enough  to 
prevent  the  overturning  tendency  induced  by  the  shear  flow. 

If  we  take  the  limit  as  R  the  eigenvalue 

equation  for  C  takes  the  form: 

[(U,  -  C)^k  -  g](AC*k  +  (U,  -  C)^k  -  g']  *  0  (23) 


The  first  factor  represents  surface  gravity  waves  being 
advected  by  the  mean  flow  (C  s  U|  ±  (g/k)^) ,  while  the  second 
factor  represents  the  Kelvin-Helmholtz  result: 


C  = 


Pi^i  ^  ^  g  Pz  -  Pi 

p,  +  p2  k  P2  +  p. 


PiPz 


(P,  Pz)' 


(24> 


or 
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2  ^1^2  2  1/2 

C  =  U  ±  (C/ - - u,  )  =  U  ±  s  (25) 

(p,  +  Pg)* 

This  result  shows  interfacial  short  waves  moving  relative  to 
the  weighted  mean  flow  U  with  speed  ±8,  provided  that  s^>0;  if 
not,  C  is  complex  and  the  flow  becomes  unstable.  This 
condition  is  valid  for  small  density  differences  and  the 
stability  criteria  is: 


Rk-  =  g^ll^  >  I 

PiU,"  2 


(26) 


(and  <  1/2  for  instability) .  The  value  Rk*'  represents  a 
Richardson  number  (Ri) ,  based  on  the  wavenumber  of  the 
disturbance.  This  inequality  reflects  the  stabilizing  effect 
of  stratification  that  is  required  in  order  to  inhibit  the 
overturning  tendency  introduced  by  the  vertical  shear. 

Now  consider  a  continuously  stratified  fluid  and 
take  the  Taylor-Goldstein  equation  assuming  a  rigid  lid  at  the 
surface,  M*=gpVP  (Brunt-Vaisaala  frequency),  w=FW,  IfcU-C, 
where  F(z)  is  a  new  dependent  variable,  and  the  following 
equation  is  obtained: 

(PqW^f')'  +  po(N^  -  kV)F  =  0  (27) 
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Now  if  f  is  an  unstable  solution,  C  is  complex  and  W  #  0. 
Using  the  substitution  FsO/W^^,  where  O  is  another  independent 
variable  and  depth  integrating  the  equation,  the  following 
condition  for  a  zero  imaginary  part  is  obtained: 

C,[JPo(|g'|*  Jpol~  I  +  N^)]  =0  (28) 


Therefore,  C|  >  0  (unstable  waves)  implies  <  0  for 

some  range  of  s  i.e.,  for  0*^0  the  sufficient  condition  for 
stzdsility  on  a  stratified  shear  flow  is: 

1 

R,=— >-  (28) 

U*  4 

known  as  the  Miles'  theorem  for  stability.  It  states  that  if 
the  stabilizing  influence  of  stratification  dominates  the 
destabil  ing  influence  of  the  shear,  the  flow  is  stable. 

From  the  linearized  perturbation  equations,  the 
following  energy  budget  is  obtained  (p.  409) : 

— {-Po(u^  +  w^)  +  -PqN^C^  =  -(PU);,  -  (pw),  -  PqUwu'  (30) 
Dt  2  2 

where  u  is  the  horizontal  velocity  and  p  the  pressure  and  C  is 
the  vertical  particle  displacement.  Assuming  that  the 
quantities  are  periodic  in  the  horizontal,  and  averaging  over 
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one  wavelength,  the  depth  integrated  energy  transfer  equation 
is  obtained  (horizontal  averages  will  be  denoted  by  B[]): 


d 

—  (T  ♦  V)  »  P  +  Q 

at 


where 


T  = 


+  w*]dz 


is  the  perturbation  kinetic  energy. 


'2 


is  the  perturbation  potential  energy 

^2 

P  =  -pwj 


represents  the  work  done  by  the  external  pressure 
boundaries ,  and 

Q  = 


(31) 


(32) 


(33) 


(34) 


at  the 


(35) 
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is  the  rate  of  energy  transferred  between  the  mean  and  the 
perturbed  flow  by  the  Reynolds  stress: 

T„  »  -PoE[uw]  (3C) 

This  equation  shows  that  when  an  instability  occurs, 
([T4>V],>0},  energy  is  transferred  from  the  mean  shear  towards 
the  perturbation  flow,  resulting  in  a  loss  of  mean  kinetic 
energy. 

In  fact.  Holt,  et  al.  (1992)  using  direct  numerical 
simulations  observed  that,  for  a  Richardson  number  (Ri)  below 
critical,  the  Reynolds  stress  and  vertical  density  flux  are 
down-gradient,  and  some  of  the  vertical  kinetic  energy  gained 
through  production  is  transferred  to  potential  energy.  As  Ri 
becomes  higher  than  critical,  there  is  a  release  of  potential 
energy  towards  the  vertical  kinetic  energy.  Some  of  this 
energy  is  then  redistributed  to  the  streamwise  kinetic  energy. 

There  are  other  mechanisms  by  which  internal  waves 
can  break  (Leblond  and  Mysak,  1978,  p.  408).  The  dynamical 
gravitational  instability,  in  which  the  vertical  density 
gradient  associated  with  large  amplitude  internal  waves 
becomes  positive,  is  a  non-linear  instability  that  occurs  when 
the  particle  velocity  exceeds  the  wave  phase  velocity.  Another 
possibility  is  the  development  of  turbulence  from  decreases  of 
the  density  gradient  due  to  local  steepening  of  isopycnals. 
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Finally,  there  are  some  experimental  results  showing 
overturning  occurring  when  large  amplitude  waves  encounter  a 
discontinuity  in  the  density  profile. 

B.  Problem 

Once  the  main  source  of  vertical  shear  is  found, 
shear  instabilities  can  be  predicted  for  real  or  idealized 
stratification  profiles  in  the  surveyed  area  using  Miles' 
criteria  (Ri  >  1/4)  for  stability,  and  compared  with  the 
occurrence  of  internally  homogenized  layers  in  our 

observations.  Also,  by  estimating  the  amplitudes  and  phase 
velocities  of  the  internal  tides,  potential  dynamical 
instabilities  can  be  identified, 
c.  Approach 

Gradient  Richardson  nuxsber  criteria  are  calculated 
for  the  instantaneous  internal  tide  baroclinic  field  that  is 
modeled  in  the  next  chapter  to  estimate  areas  where  mixing  are 
likely  to  occur.  The  number  of  occurrences  of  critical 
conditions  during  a  24  hour  cycle  are  mapped  for  a  continental 
slope,  similar  to  MBAY  conditions,  and  comparison  with  the 
observations  are  made. 
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ZV.  TKBORBTZCAL  APPROACH 


A.  ZMTERNAL  WAVE  GBMBRATZOM  (BAINES,  1982) 

This  section  describes  the  internal  tide  generation  model 
of  Baines  (1982) ,  applied  to  conditions  typical  of  the 
continental  shelf  close  to  Monterey  Bay,  to  assess  the 
importance  of  internal  tide  generation  processes  in  the  MBAY 
area,  and  to  compare  with  the  MBAY  observations.  A  simplified 
2~D  continental  slope  geometry  and  barotropic  tide  will  then 
be  used  to  estimate  tidally  forced  internal  waves  at  the  shelf 
break.  These  internal  wave  sources  will  then  be  used  to  force 
the  nonlinear  internal  wave  propagation  model  developed  in  the 
next  section. 

The  governing  equations  for  the  internal  tide  generation 
problem  are: 

au  -  - 

PoT-'^Po  fxU+V^+pgz  =  0 
ot 

dp  3po  (37) 

—  +w — -  =  0 

at  az 
v.u  =  0 

where  the  overlines  represent  vector  quantities,  'U*  is  the 
velocity  field,  'pg*  is  the  density  in  static  equilibrium  and 
■p*,  *p*  are  the  perturbed  density  and  pressure  relative  to 
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the  hydrostatic  equilibrium.  Assuming  that  the  velocity  and 
pressure  field  can  be  separated  into  barotropic  and  internal 
tide  components,  these  equations  can  be  written  in  the 
following  form: 


dU, _ 1 

— ■+  fxu.+  — Vb,+ 
at  '  Po  ' 


P/  gz  -  gp,_ 

- =  F  =  - - Z 

Po  Po 

^P,  ^Po 

It  ‘ 


(38) 


Where  *P'  represents  the  generation  force  for  internal  motion. 
The  amplitude  of  the  forcing  term  can  be  written  in  the 
following  way,  assuming  one  dominant  component  in  the  time 
dependence  of  the  barotropic  tide: 

F  .  -  9  w  ,  ^  w  ,3„ 

Po  dz  oy  <Oj. 

where  'br*  is  the  principal  tidal  frequency  that  will  be 
assumed  to  be  2.ir/12.42hr*^=1.4xlo'*8*^.  The  Brunt -Vaisala 
frequency  is  maximum  at  the  thermocline  with  magnitudes  of  the 
order  of  2xl0'*8’\  decreasing  as  the  stratification  weakens 
with  depth  towards  values  of  the  order  of  2x10  s'  for  the  MBAY 
case  (see  Table  IV-A.l).  The  semidiurnal  forcing  frequency  b^- 
is  much  smaller  than  the  Brunt-Vaisala  frequency  close  to  the 
thermocline  and  of  the  same  order  of  magnitude  as  the  inertial 
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frequency  *f' ,  which  is  O.STxlo’^s'^  for  the  latitude  of 
Monterey  Bay;  therefore,  solutions  to  this  problem  can  be 
assumed  to  be  dominated  by  near-inertial  internal  waves 
(Kunze,  1985) . 

For  a  typical  2-D  continental  slope  shallowing  eastward 
(positive  x-direction)  with  a  topographical  length  much 
smaller  than  the  barotropic  tide  scales  (-200Km) ,  the 
propagation  aspect  of  the  barotropic  tide  can  be  neglected, 
and  the  tidal  flux  can  be  assumed  to  be  constant.  The  scaling 
conditions  are  well  met  for  Monterey  Bay  even  though  the  2-D 
approximation  oversimplifies  the  topography.  By  energy 
conservation,  the  amplitude  of  the  horizontal  barotropic 
velocities  over  zhe  slope  for  a  homogeneous  ocean  for  these 
conditions  is  given  by: 

U-.  hoo 

U(x)  =  - -  (40) 

h{x) 


where  'Dj*  is  the  oceanic  tide  velocity  amplitude,  'h,'  is  the 
depth  offshore  and  'b(x)<  the  depth  over  the  continental 
slope.  Using  continuity,  with  no  lateral  boundaries  and  no 
normal  flow  perpendicular  to  the  vertical  boundaries,  the 
barotropic  velocity  interacting  with  the  continental  slope 
forces  vertical  motions  with  amplitudes  given  by: 
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(41) 


dh  z 

W(x, z)  *  U(x) - 

dx  h(x) 

Representative  values  in  the  neighborhood  of  the  surveyed  area 
during  the  MBAY  experiment  were  U-isO.OSm/s,  dh/dxso.l8  (slope) 
with  depths  ranging  from  h(0)sioom  to  b.=1000m  over  a  5000m 
horizontal  distance  across  the  slope.  Figure  IV-A.l.a/b 
represents  x/z  sections  of  the  resultant  velocities  for  the 
cross-shore  (u)  and  vertical (w)  velocity  components.  For  a 
strong  thermocline  as  observed  during  MBAY,  this  velocity 
field,  will  produce  a  maximum  internal  tide  forcing  close  to, 
but  below  the  thermocline  at  the  shelf  break,  as  Figure  IV- 
A.l.c  shows. 

Defining  a  stream  function  (f )  as  U|B-tz  where  the 

subscripts  »x»  and  represent  partial  derivatives,  a  forced 
wave  equation  can  be  derived: 

*  QN^z(-)„cos(avt)  (42) 

h 


where  *Q*  represents  the  tidal  flux  assumed  to  be  For 

a  typical  solution  of  the  form  exp(-ibjt)  we  obtain: 


(43) 
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where  the  boundary  conditions  are  T=0  at  s=0  and  -h(z),  and 
*o'  represents  the  semidiurnal  internal  wave  slope  (ratio  of 
horizontal  and  vertical  wavenumbers)  given  by: 


c 


2 


(44) 


The  internal  wave  slope  for  the  MBAY  data  will  be  minimum  at 
the  thermocline  with  a  magnitude  on  the  order  of  0.006  for  an 
interface  temperature  drop  of  3®c  (Figure  III-B.l)  ,  increasing 
with  depth  to  0.06  as  'M*  approaches  2»10'^a*\  below  the 
thermocline  (see  Table  IV-A.l). 

A  soli  tion  for  the  homogeneous  problem  can  be  generated  in 
terms  of  an  infinite  set  of  normal  modes,  assuming  a  density 
interface  at  depth  '-d',  with  a  constant  'Ng*  below,  and  a 
constant  depth.  These  modes  will  be  determined  by  the 
eigenvalues  satisfying  the  relation: 

(K^^S  -  K)tan(K„(l  -  R))  =  K„R  (45) 


where 
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(46) 


R 

S 


d/h 

gApo 


=  0.2 

/  (No  -  4) 


dPc 


K  = 


K„RC/C, 


tanh(K„Rc/c,) 

1 

f  2.  2 


C,  -  (1  -  (  — )") 


=  35 


These  values  are  representative  of  the  MBAY  conditions  at  the 
shelf  break  with  ds20n  and  hslOOm,  and  where  Ng  was  estimated 
based  on  the  temperature  maps  shown  in  Figure  III-B.l  as 
2.9*10"*s’\  corresponding  to  a  0.25®C  temperature  drop  over  a 
50m  depth,  and  the  temperature  difference  across  the 
thermocline  interface  is  estimated  to  be  3*C.  For  these 
values: 


R  > 


0.097 


(47) 
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implying  that  for  the  NBAY  conditions  only  the  lower  mode  will 
be  significant  at  the  interface  depth  (Baines,  1982). 
Therefore,  an  interfacial  mode  can  be  expected,  decoupling  the 
upper  layer  from  the  internal  wave  propagation  in  the 
stratified  continuum  below.  This  result  allows  the  two  layers 
to  be  treated  separately  which  is  consistent  with  the 
observations  discussed  in  III-B.  The  observations  showed 
different  responses  in  the  layer  above  the  main  thermocline 
compared  to  the  layer  below,  although  both  were  apparently 
tidally  forced. 

For  the  case  of  the  interfacial  mode,  the  above  wave 
equation  can  be  integrated  across  the  density  interface  to 
obtain: 


Qd 

T(x,  -  d)  »  -  —  +  0^ 
h  ^ 


(48) 


and 


-  ox  - 

d 

h  +  z 

♦^(x,  -  d)  ( - ) 

^  h  -  d 


-  d  <  z  <  0 

-  h  <  z  <  -  d 


(49) 


where 
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X 


X  >  0 


(50) 


♦  =  R  +  D^e 


(1  -fl) 


2Hq  2Mo 

+  (1  -  — T)#  - - 

r  r 


X  <  0 


where  Dg  is  obtained  by  the  boundai  .ions  (X=0)  at  the 

shelf  break  and 


3000ia 


Q 

Xg  =  T(1  -  R)  ; 
r  =  Xg  +  X 


(51) 


and  g*  represents  the  reduced  gravity,  h|^  is  the  depth  at  the 
shelf  break  (lOOm  for  Monterey  Bay),  a  is  the  bottom  slope 
(0.18  for  the  simplified  2-D  geometry  being  used)  and  2irL  is 
the  wavelength  of  the  inter facial  mode  (app.  18000m  for  the 
MBAY  case)  .  Typical  solutions  for  the  equation  at  X  <  0  (i.e., 
over  the  slope)  can  be  found  in  terms  of  expansions  of  the 
Coulomb  wave  functions  (Abramowitz,  Stegun,  1968) . 

In  order  to  estimate  solutions  over  a  steep  slope,  a  step 
geometry  is  considered  using  typical  MBAY  parameter  values  as 
shown  in  Table  III-A.l  and  corrections  for  the  slope  case 
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determined.  Therefore,  for  the  step  geometry  case,  the 
following  results  are  obtained: 


« 


=  R  +  D^e" '  ‘  ^ 


#X(1  -%) 


J 

1 

1 
’  -i 


X  >  0 

X  <  0 


(52) 


Where  R||=d/h«.  From  the  boundary  conditions  we  obtain: 


R  -  Rb 


Do  = 


1  -  Ro 

1  ^ 
1  -  R 


D,  =  -  ( - -)  \ 

1  -  R 


(53) 


For  the  MBAY  conditions  OoSO.085  and  D)S0.095;  however,  Baines 
(1982)  shows  that  corrections  for  a  linear  continental  slope, 

"I0*2v  "Id 

changes  these  constants  to  be  DgsO.oese  and  D^so.076e^, 
where  'p*  is  a  phase  factor  in  the  Coulomb  wave  functions  and 
is  determined  by  depth,  slope  and  *T*.  Furthermore,  using 
expressions  (7.4)  from  Baines  (1982)  the  current  magnitude 
over  the  shelf  is  estimated  to  be  0.17m/s  and  0.04m/s  in  the 
upper  and  lower  layers  respectively,  with  an  interfacial 
displacement  amplitude  of  8.9m.  These  values  will  not  change 
significantly  over  the  continental  slope  since  T  is  small, 
resulting  in  0.19m/s  and  0.004m/s  for  the  slope  current 
magnitudes  in  the  upper  and  lower  layer  respectively  and 
interface  displacements  of  8.9m  (similar  to  the  shelf  values)  . 
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However,  the  phase  lag  is  not  easily  determined  because  of  the 
phase  variation  of  the  Coulomb  functions  with  R,  T  and 
alongslope  distance.  As  a  result,  for  the  MBAY  case  siunmarized 
in  Table  IV-A.l,  a  maximum  velocity  gradient  of  0.20m/s  can  be 
expected  across  the  thermocline  due  to  the  interfacial  mode, 
resulting  in  a  Richardson  number  of  about  1.4  for  a  10m  width 
thermocline,  i.e.,  small  but  above  the  critical  value  of  0.25. 
Therefore,  the  thermocline  oscillations  with  ~10m  amplitudes, 
wavelengths  larger  than  lOOOOm,  and  tidal  periodicity  which 
were  observed  during  the  first  part  of  the  survey  and 
discussed  in  Chapter  III-B,  are  predicted  using  this  theory. 
Furthermore,  the  wavelength  of  this  interfacial  mode  will  be 
of  the  order  of  18000m,  with  an  e-folding  distance  of  35  Km; 
consequently  this  mode  will  have  significant  amplitudes  at 
considerable  distances  from  the  shelf  break.  This  theory 
partly  explains  the  different  observed  velocity  fields  in  the 
layers  above  and  below  the  thermocline  until  hour  152, 
discussed  in  Chapter  III-B.  However,  the  velocity  fields  in 
both  the  upper  and  lower  layers  displayed  in  Figures  III-B.  4/5 
at  11m  and  50m  depths  do  not  show  a  clear  tidal  signature, 
implying  that  other  baroclinic  phenomena  are  also  present. 

From  the  results  of  Baines  (1973)  ,  the  internal  wave 
amplitudes  in  the  stratified  continuum  below  will  be  small 
unless: 
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(54) 


1  -  a/c 

p  »  - 1_  <  — 

1  +  o/c  h« 

This  inequality  is  observed  in  the  MBAY  case  so  significant 
amplitude  internal  waves  can  be  expected  to  be  generated  below 
the  thermocline.  In  fact,  using  expression  (7.9)  from  Baines, 
(1982) ,  internal  wave  amplitudes  can  be  estimated  by: 


A  »  - 


h  -  d 


(1  -  D,)  - 


4  h^^  -  d 


(1  -  D,) 


“r-  ^  -2 


- 


X[V(  — )  -V(— 1)] 
2y  2y 


(55) 


where  V<{)  is  the  function  developed  in  Baines  (1982) 
characterizing  the  wave  amplitudes  across  the  internal  wave 
beam  generated  at  the  shelf  break.  The  variables  i,  rj  and  y 
are  determined  by  expressions  (7.10)  from  Baines  (1982) 
representing  the  normalized  distance  across  the  beam,  distance 
along  the  beam,  and  wavelength  scale  respectively.  Therefore, 
assuming  a  uniform  stratification  below  the  thermocline  with 
Ms2.9xl0  8*  for  a  typical  characteristic  inside  the  beam, 
V({/2y)"V(-»)/2y)*”2  f  results  in  a  baroclinic  velocity 
amplitude  O(0.15m/8).  Assuming  with  the  same  phase  as  the 
step  geometry  solution,  these  oscillations  will  lag  the 
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barotropic  tide  by  9/2,  and  will  have  a  typical  wavelength 
given  by 


.  ^ 

7 1 

Y  =  — dz  =  2000  m 


(56) 


However,  because  the  stratification  is  not  uniform  but 
concentrated  in  a  shallow  thermocline  in  the  MBAY 
observations,  the  resulting  forcing  is  dominated  by  a  shallow 
layer  below  the  thermocline  at  the  shelf  break. 

B.  ZMTBBIIAL  WAVE  PKOPAGATIOM 

This  section  summarizes  the  numerical  procedure  used  to 
study  the  non-linear  propagation  of  the  internal  waves 
generated  at  the  shelf  break  interacting  with  the  barotropic 
tide.  For  a  better  understanding  and  interpretation  of  the 
results,  a  single  component  wave  solution  will  be  propagated 
from  a  single  point  source  ar  the  shelf  break,  where  the 
forcing  is  predicted  to  be  a  maximum.  Therefore,  the  following 
results  can  be  interpreted  as  an  elemental  response  that  might 
be  integrated  later  with  the  solutions  starting  from  other 
points  within  the  water  column,  so  that  a  more  general 
solution  to  the  problem  can  be  obtained.  The  focus  of  this 
method  will  be  to  simulate  the  nonlinear  interactions  of 
internal  waves  induced  by  the  barotropic  tidal  velocities, 
which  are  modified  by  topography  along  a  2-dimensional 
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continental  slope.  This  approach  will  investigate  the 
potential  of  these  phenomena  to  generate  the  strong, 
baroclinic,  jet-like  pulses  observed  during  the  NBAY 
experiment  and  discussed  in  Chapter  III-C. 

Assuming  a  mean  flow  field  varying  slowly  in  time  and 
space,  within  which  a  narrow  bandwidth  internal  wave  is 
propagating,  the  Boussinesq  approximation  to  the  equations  of 
motion  (Leblond,  Mysak,  1978)  can  be  written  in  the  following 
form,  neglecting  wave-wave  interaction  terms  (Kunze,  1985) : 

u^  +  (V.V)u  +  (v.V)U  -  f.v  =  -  Pj, 

Vj  +  (V.V)v  +  (v.V)V  -t-f  .u  =  -  Py 

Wj  +  (V.V)w  +  (v.V)W  =  -  p^  -  b  <*7) 

u,  Vy  w,  =  0 

b/  +  (V.V)b  +  uB^  +  vBy  -  *  0 

where  the  (u,v,w)  represent  the  internal  wave  velocities  and 
(U,V,W)  the  mean  flow  velocities.  Terms  (V.V)v  correspond  to 
the  advection  of  the  wave  solutions  by  the  underlying  flow  and 
are  hereafter  called  'advective*  or  'Doppler  shifting*  terms. 
The  terms  (v.V)y  correspond  to  the  straining  of  solutions  by 
the  underlying  flow  and  are  hereafter  called  'straining' 
terms.  The  mean  flow  and  wave  buoyancy  are  represented  by  'B* 
and  *b'  respectively,  and  *p*  is  the  wave  induced  pressure 
normalized  by  the  mean  density  of  the  Boussinesq  fluid.  The 
mean  flow  is  assumed  to  satisfy  the  following  relations: 
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(58) 


U,  -  f .V  =  - 

V,  +  f  .U  -  - 

W,  =  -  P,  -  B 


SO  that 


Bx  =  -f-V,  -  (U,  -W^), 
By  =  f-U^  -  (V,  -  W^), 


(59) 


If  we  assiune  that  an  elemental  internal  wave  variable  can 
be  written  in  the  form  of  a  plane  wave 

T  =  i;,(x,t)et''<'^‘’'^^  (50) 

where  the  rate  of  change  of  amplitude  is  small  compared  to  the 
phase  variations,  then  the  equations  of  motion  can  be  written 
in  terms  of  the  wavenumbers  K  and  of  the  intrinsic  frequency 
following  the  mean  motion  ttgSe-K.V,  where  s  is  the  Eulerian 


frequency  and  K.V  the  doppler  shifting: 


-ieo+Dx 

-f  +  Uy 

Bx 

i  Kx 

0 

u 

= 

0 

f  +  V, 

-i«0+Vy 

▼x 

i  Ky 

0 

V 

0 

1  Xs 

1 

w 

0 

-'▼x+C/ 

0 

-itto 

p 

0 

i  Xz 

i  Ky 

i  Ks 

0 

0 

b 

0 

69 


with,  C^sVj-Wy  and  C"*Wy-0„  and  (Kx,Ky,K8)  represents  the 
wavenuaber  components.  To  obtain  non-trivial  solutions,  the 
determinant  of  this  system  must  be  zero.  The  roots  of  the 
resultant  3'*^  order  polynomial  with  complex  coefficients  form 
an  implicit  dispersion  relation  which  will  be  used  to 
numerically  trace  the  ray  path  (characteristic  curve  of  the 
system)  and  predict  the  frequency  and  wavenumber  of  an 
elemental  solution.  For  a  3-D  case,  no  explicit  polynomial  was 
obtained  and  the  roots  were  computed  numerically  from  the 
implicit  relation. 

Typically,  the  roots  of  this  characteristic  polynomial 
represent  three  modes:  one  purely  complex  corresponding  to  the 
vortical  mode;  and  the  other  two  with  the  same  imaginary  part 
but  having  symmetric  real  parts  representing  the  internal 
gravity  modes  (Lamb  et  al.,  1992).  However  other  solutions  are 
possible  depending  on  the  wave  group  and  phase  speeds  relative 
to  the  mean  flow.  For  instance,  for  evanescent  waves  all  the 
roots  are  imaginary  (Lamb  et  al.,  1992),  and  when  the  waves 
and  mean  motion  are  in  opposite  directions,  non-symmetric 
roots  with  only  negative  real  parts  are  possible,  resulting  in 
the  waves  being  advected  faster  then  their  phase  propagation 
as  mentioned  by  Peregrine  (1975)  for  surface  gravity  waves. 

The  solutions  to  this  dispersion  relation  will  then  be 
used  following  an  approach  usually  referred  as  the  WKB 
approximation  (Leblond  and  Mysak,  1978;  Bretherton  and  Garret, 
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1968;  Bretherton  1971;  Bender  and  Orzag,  1978),  which  has  been 
used  by  Peregrine  (1975)  to  study  surface  gravity  waves 
interacting  with  the  mean  flow.  Baines  (1982)  ,  as  mentioned  in 
the  previous  section,  used  an  equivalent  approach  to  study  the 
propagation  of  internal  tides  in  continuously  stratified 
fluids;  however,  he  did  not  consider  the  non-linear  effects  of 
the  barotropic  tide  (doppler  shifting  and  straining) .  Using 
the  same  linear  approach,  Baines  and  Fang  (1985)  compared  the 
results  of  the  model  with  a  laboratory  experiment,  and 
concluded  that  non-linear  effects  should  be  important  in  a 
layered  fluid  at  high  (above  unity)  Froude  number,  Fr»»U/Cp, 
where  U  represents  the  tidal  speed  at  the  shelf  break  and  Cp 
the  phase  speed  of  the  lowest  mode.  The  MBAY  data  showed 
strong  currents  (unidirectional,  baroclinic,  alongslope  pulses 
eveiry  12-24  hours) ,  a  strong  internal  mixed  layer  below  the 
thermocline,  and  stratification  with  a  sharp  thermocline  at 
about  20m  depth,  suggesting  that  non-linear  dynamics  are 
relevant.  Furthermore,  from  the  buoy  data  of  Pillsbury  (1992) , 
barotropic  tidal  velocities  are  estimated  to  be  below  O.lOm/s. 
Based  on  the  results  of  the  previous  section,  the  phase  speed 
of  the  interfacial  mode  is  estimated  to  be  about  0.40m/s, 
resulting  in  a  Froude  number  near  the  thermocline  well  above 
unity  (-4) ,  at  least  for  some  portion  of  the  tidal  cycles. 

Willmott  (1987)  and  Maze  (1987)  developed  three  and  two 
layer  nonlinear  models  respectively,  taking  into  account 
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advection  by  the  barotropic  tides.  They  showed  that  velocities 
of  the  resultant  internal  wave  field  can  be  several  times 
higher  than  the  barotropic  tide,  and  that  the  resultant 
baroclinic  residual  flows  change  in  space  with  the  same  length 
scales  as  the  internal  waves  as  a  result  of  the  "bunching" 
effect  that  sometimes  can  halt  the  internal  wave  propagation. 
However,  they  did  not  take  into  account  rotation  and  straining 
effects  by  the  tidal  flow.  Furthermore,  they  used  two  and 
three  layered  models,  as  approximations  to  the  physical 
system,  which  do  not  necessarily  well  represent  the  baroclinic 
propagation  of  the  internal  tides. 

Based  on  the  linear  modal  model  by  Prinsenberg  and  Rattray 
(1975),  the  internal  waves  should  propagate  as  beam-like 
structures  away  from  the  shelf  break.  New  (1988)  found  tidal 
currents  along  these  beam-like  current  structures  (as 
described  also  by  Baines,  1982)  with  velocities  0(0.50m/s)  in 
excess  of  the  surface  tides  in  the  Biscay  Bay.  He  also  showed 
that  the  strong  shear  created  by  the  internal  tides  was  likely 
to  generate  strong  mixing  layers.  Kunze  (1985)  used  a  ray 
tracing  (WKB)  method  to  study  near-inertial  internal  waves 
propagating  in  geostrophic  shear,  and  showed  the  importance  of 
including  the  straining  terms  in  the  nonlinear  equations  of 
motion.  He  also  showed  that  the  WKB  method  was  still 
approximately  valid  even  when  the  spatial  scales  of  the  mean 
flow  were  comparable  to  the  wavelengths  of  the  internal  waves. 
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The  present  developnent  extends  Kunze's  work  by  considering 
the  time  variability  in  the  underlying  mean  flow. 

The  ray  equations  for  wave  propagation  in  a  moving  media, 
using  indicial  notation,  (Leblond  and  Mysak,  1978)  are: 

d  X, 

—  *  —  +  u, 

d  t  dK, 

d  K/  _ 

d  t  dX  dXj  ‘  dXf 

where  the  total  derivatives  represent  rates  of  change 
following  the  ray  path  and 

<0  =  o^,[K(x,t)  ;X(x,t)  ]  +  K.U  («) 


is  the  dispersion  relation  resulting  from  the  equations  of 
motion  discussed  above,  with  A(x,t)  representing  properties  of 
the  medium  stratification  and  straining  by  mean  flow. 
Frequency  is  not  constant  along  characteristics  due  to  the 
time  variability  of  the  mean  flow,  and  an  explicit  form  of  the 
dispersion  relation  is  required  to  compute  the  derivatives  in 
the  ray  equations.  In  the  present  work,  it  was  not  possible  to 
obtain  an  explicit  relation  for  a  3-D  case,  so  the  implicit 
expression  represented  by  the  roots  of  the  3'*^  order  polynomial 
described  above  are  used  to  numerically  compute  the  required 
derivatives.  The  largest  real  part  of  the  roots  was  chosen  as 
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the  physical  solution  to  the  dispersion  relation.  When  all  the 
roots  are  imaginary  or  have  negative  real  parts,  the  wave 
solution  is  assumed  evanescent  with  all  its  energy  absorbed  in 
some  way. 

The  wave  amplitudes  following  characteristics  starting  at 
the  shelf  break  and  travelling  along-slope  can  be  estimated  by 
considering  the  conservation  of  wave  action  'A* ,  defined  as 
the  ratio  of  the  mean  wave  energy  to  the  intrinsic  frequency 
as  follows  (Bretherton,  1968;  Garret,  1968;  Bretherton, 
Garret,  1969;  Peregrine,  1975;  Leblond,  Mysak  1978): 


dA 

dt 


-V.C,  A 


«3) 


where  'C^'  is  the  group  velocity  defined  as 

*  — 

®  dKi 


(64) 


The  divergence  of  the  group  velocity  is  computed  including  the 
effects  of  spatial  inhomogeneities  of  the  mean  flow  and 
stratification,  and  the  effects  of  spatial  changes  of  the 
wavenumbers.  In  indicial  notation  this  becomes: 


v.c. 


_ £' _ >  I  ^ 

3k.  dt  C  .  +  u.  *"®®"**’  3X. 


K^onst. ' 


i  =  1,  2  (65) 
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The  implementation  of  this  approach  is  summarized  next 
using  the  linear  continental  slope  and  simplified  tidal  flow 
described  in  the  previous  section,  and  assuming  an  internal 
wave  solution  at  the  shelf  break  as  predicted  in  the  previous 
section  using  the  theory  of  Baines  (1982) .  A  plane  wave  is 
assumed  as  a  typical  solution  to  a  first  order  approximation, 
allowing  the  characteristic  paths  to  be  estimated  using  a  ray 
tracing  method.  A  known  initial  solution  is  taken  at  the  shelf 
break  (x=0)  where,  based  on  the  results  of  the  last  section, 
a  single  resonant  Fourier  component  type  of  oscillation  can  be 
assumed  to  exist  (Haberman,  1987;  Leblond  et  al.,  1978; 
Witham,  1974)  .  This  initial  solution  is  idealized  to  have  non¬ 
zero  amplitude  only  at  a  single  depth  close  to  the  thermocline 
(-30m  depth)  and  over  the  shelf  break  (-zero  width)  and 
follows  the  shape  of  the  forcing  function  displayed  in  Figure 
IV-A.l.c.  Therefore,  this  analysis  considers  characteristics 
originating  only  from  a  single  point  in  the  water  column, 
located  at  the  shelf  break  and  at  30m  depth.  However,  because 
the  initial  solution  and  the  environmental  conditions  change 
in  time  following  the  tidal  cycle  (which  represents  the  non¬ 
stationary  mean  flow) ,  the  characteristics  at  each  time  in  the 
tidal  cycle  need  to  be  computed  at  that  origin  (Witham,  1974) . 
The  solution  at  the  origin  will  propagate  along  the 
characteristic  lines,  so  that  each  point  along  the  line  will 
have  a  time  lag  relative  to  the  initial  condition,  and  a 
wavenumber  determined  by  the  changes  in  the  mean  flow.  The 
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wavelengths  increase/dec  'ease  as  the  wave  propagates  into  a 
faster/slower  flow  (Peregrine,  1975) .  The  next  step  uses  the 
wave  action  conservation  (Peregrine,  1975  p.37)  to  compute 
velocity  changes  along  these  paths.  These  estimates,  together 
with  the  time  lag  and  initial  phase,  allow  the  computation  of 
the  instantaneous  velocities  of  the  solution  at  each  point 
along  the  ray  path.  The  velocity  estimates  correspond  to  the 
square  root  of  half  the  wave  action  times  the  ray  Euler ian 
frequency.  The  velocity  solution  changes  in  time,  following 
the  instantaneous  Eulerian  frequency  estimates,  which 
increase/decrease  as  the  flow  accelerates/decelerates 
(Peregrine,  1975) .  These  changes  in  magnitude  and  frequency 
make  the  solutions  nonlinear,  inducing  different  residual 
flows  along  the  characteristics. 

It  is  assumed  that  reflections  at  the  surface  and  bottom 
have  no  losses  and  that  there  is  no  dissipation  by  frictional 
effects.  These  assximptions  could  be  easily  changed  in  the 
future  to  simulate  more  realistic  environments.  Rays  with  the 
intrinsic  frequency  reaching  zero  or  with  the  wavenumbers  and 
wave-action  increasing  toward  non-physical  values  are  assumed 
to  be  dissipated  by  some  type  of  instability,  such  as  critical 
layer  processes,  other  areas  where  the  rays  are  assumed  to  be 
dissipated  are  caustics  i.e.,  rays  converging  toward  a  line  or 
point  in  space  with  no  possibility  of  reflection,  where  the 
WKB  or  the  plane  wave  approximation  is  completely  invalid 
(Peregrine,  1975) .  The  ray  equations  and  the  wave  action 
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balance  are  integrated  using  a  4  order  Runge-Kutta  scheme  with 
an  adaptative  time  step  satisfying  a  specified  threshold 
criteria.  The  time  steps  were  all  less  than  1  hour,  and  close 
to  surface  or  bottom  reflections  were  a  few  seconds. 

To  verify  the  numerical  procedure  independently  of  the 
implicit  limitations  of  the  plane  wave  assumption  that  will  be 
discussed  below,  a  2 -dimensional  formulation  with  explicit 
expressions  for  the  dispersion  relation  is  also  used.  The 
results  for  the  steady-state  mean  flow  developed  by  Kunze 
(1985)  are  also  used  for  a  qualitative  comparison  of  the  3-D 
model  runs.  All  the  numerically  computed  terms  in  the  ray 
equations  are  compared  with  the  analytical  2-D  counterparts 
for  a  typical  ray.  They  showed  small  departures  on  the  order 
of  10%  and  similar  characteristics  for  the  first  30  hours, 
except  for  the  group  velocity  divergence  which  show  larger 
differences.  These  differences  are  dominated  by  the 
introduction  of  the  coriolis  term  in  the  dispersion  relation. 
However,  as  the  group  velocity  divergence  requires  the 
numerical  computation  of  second  derivatives,  which  amplifies 
numerical  noise,  it  was  decided  to  introduce  a  wave  action 
maximum  threshold  value  of  30  times  greater  than  the  initial 
value,  after  which  characteristics  exceeding  this  threshold 
were  terminated  (physically,  the  ray  was  dissipated) .  This 
limit  threshold  value  was  chosen  based  on  the  wave  action 
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estimates  for  typical  straight  downward  rays,  which  are  in  the 
range  0.1  to  1.5,  times  the  initial  wave  action  value. 

To  validate  this  type  of  approach  for  the  case  of 
propagation  of  near-inertial  internal  waves  on  topographically 
influenced  tidal  flow,  several  points  must  be  considered. 
Usually  the  mean  flow  is  assumed  to  be  slowly  changing  and  to 
have  length  scales  much  larger  than  the  internal  wave  length 
scales,  so  that  a  formal  expansion  in  terms  of  the  following 
small  parameters  can  be  made  (Bretherthon,  1968;  Bretherthon, 
1969;  Peregrine,  1975); 


,  1  dU, 

MAX - 

U8x 

€  =  -  <  1 


,  1  dU, 

MAX - 

udt 

p  =  -  «  1 


(«6) 


However,  for  the  MBAY  case  and  c»l  over  the  continental 
slope,  suggesting  that  this  approach  is  not  valid,  and  higher 
order  corrections  are  needed.  One  of  the  direct  results  from 
incorporating  these  higher  order  corrections  is  the 
possibility  of  adding  the  residual  currents  resulting  from 
previous  internal  wave  interactions  into  the  underlying  mean 
flow,  as  discussed  by  Bretherthon  (1969)  and  Peregrine  (1975) . 
However,  the  results  of  Bender  and  Orzag  (1978)  showed 
surprisingly  accurate  results  comparing  the  WKB  approach  with 
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analytical  solutions  for  cases  close  to  the  WKB  limit.  Kunze 
(1985)  also  obtained  accurate,  physically  well-based  results, 
and  successfully  compared  the  effects  of  a  critical  layer 
computed  numerically  by  a  different  method  with  the  WKB 
approximation,  demonstrating  the  presence  of  a  caustic. 

The  present  formulation  allows  time  variability  for  the 
mean  flow,  so  that  additional  factors  and  limitations  should 
be  considered.  Bretherthon  and  Garret  (1969)  suggested  that 
terms  contributing  to  the  dispersion  relation  must  be  slowly 
varying.  For  internal  waves,  they  also  suggested  a  relevant 
condition  is  that  the  ratio  of  the  horizontal  to  the  vertical 
length  scales  of  the  underlying  flow  should  be  much  larger 
than  the  ratio  of  the  internal  wave  horizontal  to  a  vertical 
wavelength.  Considering  a  barotropic  underlying  flow  with 
horizontal  length  scales  on  the  order  of  200Km,  typical  of  the 
MBAY  site,  this  last  condition  is  satisfied.  These 
considerations  motivated  the  current  approach;  however,  as 
comparison  of  these  model  results  with  other  numerical  or 
theoretical  solutions  will  not  be  possible,  validation  of  this 
procedure  is  based  on  comparison  with  observational  data  and 
physical  interpretations. 

C.  NUMERZCAL  RESULTS 

This  section  presents  the  results  of  the  implementation 
of  the  niunerical  procedure  developed  in  the  previous  section 
for  a  linear  continental  slope,  simplified  barotropic  tidal 
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flow,  and  assuming  internal  waves  with  amplitudes  and 
wavelengths  predicted  above  using  the  theory  of  Baines (1982) 
at  the  shelf  break. 

Discussion  of  the  results  of  two  different  numerical  model 
test  runs  follows.  In  run  #1,  the  numerical  model  is  forced 
with  internal  waves  with  a  frequency  six  times  the  tidal 
frequency,  horizontal  wavelengths  of  2Km  (small  compared  to 
the  topographical  scale) ,  and  a  reduced  continental  slope 
extending  over  lOKm  (compared  with  the  5Km  continental  slope 
near  Monterey  Bay) ,  and  with  depth  changing  linearly  from  200 
to  1000m.  For  this  run  the  WKB  approach  is  theoretically 
valid,  allowing  the  numerical  procedure  to  be  validated  by 
checking  the  consistency  of  the  results.  Run  #2  more  closely 
matches  the  simplified  geometry  for  the  MBAY  site  with  a 
continental  slope  extending  over  5Km  and  depth  changing 
linearly  from  200  to  1000m,  but  it  is  forced  with  internal 
waves  with  twice  the  tidal  frequency  having  a  horizontal 
wavelength  of  approximately  7Km.  This  is  similar  to  the 
topographical  scale  and  to  the  predicted  wavelength  for  the 
internal  tides,  so  the  spatial  conditions  become  nearly 
critical,  although  there  is  still  a  difference  in  time  scales. 
The  results  of  these  two  runs  are  compared  by  examining  their 
ray  paths  and  the  distribution  of  energy  along  them.  In  each 
case,  the  families  of  characteristics  are  computed  at 
different  initial  times  throughout  a  single  tidal  cycle, 
considering  both  initially  downward  and  upward  wave  packets. 
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The  reference  time  in  these  runs  (hour  0)  corresponds  to  low 
tide. 

The  analysis  of  the  model  results  compares  these  test  runs 
with  run  #3,  which  considers  the  same  geometry  as  run  #2  but 
is  forced  with  semidiurnal  internal  waves  having  7Km  length 
scales  and  amplitudes  of  0.15m/s,  as  predicted  using  the 
theory  of  Baines  (1982)  in  Chapter  IV-A.  Selected  model 
velocity  time  series  at  points  having  the  strongest  residual 
flows  are  compared  with  time  series  velocity  data  from  the 
MBAY  experiment  to  assess  the  ability  of  this  model  to 
generate  the  observed  residual  flows.  Finally,  the  number  of 
occurrences  of  Richardson  number  below  critical  (<0.25)  is 
evaluated  during  a  24  hour  period  and  mapped  to  identify  areas 
where  mixing  is  likely  to  occur.  The  results  are  discussed  and 
compared  with  the  internal  mixed  layer  observed  during  the 
MBAY  experiment.  All  these  model  runs  consider  horizontally 
homogeneous  stratification,  with  a  sharp  thermocline  at  20m  as 
summarized  in  Table  IV-C.l. 

For  the  conditions  of  run  #1,  linear  theory  predicts  that 
do%imward  rays  should  reflect  from  the  slope  approximately  IKm 
from  the  shelf  break  for  a  linear  wave  slope  of  approximately 
0.03;  however,  as  Figure  IV-C.l  shows,  the  initially  do%mward 
rays  starting  at  hours  2  and  4  reach  the  floor  closer  to  the 
shelf  break  at  distances  on  the  order  of  500m.  This  is  due  to 
the  onshore  advection  by  the  horizontal  tidal  flow  against  the 
horizontal  ray  propagation,  and  to  an  increased  wave  frequency 
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due  to  doppler  shifting,  which  further  reduces  the  internal 
wave  slope  'C  as  defined  in  Chapter  IV-A.  On  the  other  hand, 
the  straining  by  the  tidal  flow  decreases  the  wave  intrinsic 
frequency  decreasing  the  wave  momentum,  resulting  in  the 
trapping  of  these  rays  close  to  the  generation  site.  This  can 
be  seen  in  the  initial  amplitudes  plot  of  Figure  IV-C.2 
comparing  values  to  a  0.15m/s  amplitude  at  the  origin.  The 
contribution  of  these  characteristics  is  minimal  and  their 
energy  can  be  assumed  to  be  rapidly  absorbed  close  to  the 
shelf  break. 

The  downward  rays  starting  at  hour  0  and  6  experience  low 
tidal  velocity  at  the  source,  so  they  hit  the  bottom  at  a 
distance  comparable  to  that  predicted  from  linear  theory. 
However,  as  they  propagate,  they  interact  with  offshore 
velocity  fields  allowing  the  ray  starting  at  hour  6  to 
propagate  further  away.  The  intrinsic  frequency  of  these  rays 
is  a  maximum  at  the  source,  implying  by  wave  action 
conservation  that  their  momentum  is  also  largest  at  the 
source.  Therefore,  their  momentum  will  decrease  as  they 
propagate  away  from  the  shelf  break.  This  effect  appears  to  be 
stronger  for  the  ray  starting  at  hour  0,  which  is  trapped  near 
the  surface  at  3Km  from  the  source. 

The  rays  generated  at  hours  8  and  10  also  show  the  effects 
of  doppler  shifting  described  above.  The  tidal  flow  is 
decelerating  (from  onshore  to  offshore),  reducing  the  wave 
frequency  and  decreasing  the  ray  slope,  so  they  reach  the 
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floor  at  a  larger  distance  from  the  shelf  break.  The 
horizontal  advectlon  of  the  ray  at  hour  8  is  stronger  and 
directed  offshore,  so  that  the  ray  propagates  further  away. 
The  rays  starting  at  hours  0,  2  and  10,  have  accelerating  mean 
flows  at  the  origin  causing  them  to  break  down  before  others 
that  reach  the  bottom  at  the  same  distance  from  the  source. 

The  upward  rays  experience  the  same  effects  described 
above  for  waves  propagating  in  the  opposite  vertical 
direction;  however,  they  also  cross  a  layer  of  strong 
stratification  which  attenuates  their  vertical  group  velocity 
(Cg) .  Therefore,  only  rays  at  hours  8  and  10,  starting  with 
the  stronger  offshore  tidal  flow  propagate  significant 
distances  from  the  shelf  break.  For  the  present  case,  the  wave 
slope  is  larger  than  the  bottom  slope  and  the  effect  of  the 
horizontal  advection  dominates  the  ray  paths.  The  frequencies 
at  the  source  seem  to  be  important  in  the  definition  of  the 
initial  energy  distribution  and  the  distance  each  component 
propagates  from  the  origin,  so  the  rays  with  the  smaller 
initial  frequency  propagate  further  offshore. 

For  the  MBAY  case  (run  #2),  the  ray  paths  displayed  in 
Figure  IV-C.3  show  a  behavior  similar  to  run  #1,  although  the 
paths  extend  further  in  the  horizontal  as  would  be  expected 
for  this  lower  frequency.  Also,  the  increased  slope  introduces 
larger  horizontal  velocities,  further  increasing  the  range  of 
the  rays.  The  similarity  between  these  two  runs  can  be  seen  by 
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conparing  the  9Kn  distance  shotm  for  this  run  with  the  first 
2Kb  in  run  #1.  However,  one  should  note  that  more  rays 
propagate  to  greater  distances  from  the  origin,  and  a  small 
energy  enhancement  can  be  observed  in  the  downward  rays 
starting  at  hours  6-10  before  they  reach  the  bottom  (as  shown 
in  the  initial  amplitude  map  displayed  in  Figure  IV-C.4). 
Another  important  feature  is  that  the  upward  rays,  propagating 
towards  a  more  weakly  stratified  layer,  have  smaller  initial 
amplitudes  which  rapidly  change  in  space.  Some  of  the 
velocity/ amplitude  values  close  to  the  shelf  break  increase  to 
unrealistic  high  energy  levels.  This  suggests  that  the  plane 
wave  approximation  might  be  violated  for  these  cases,  and  that 
these  solutions  are  unstable.  These  results  can  be  considered 
as  a  consistent  extension  of  the  conditions  for  run  #1. 

For  run  #3,  significant  changes  are  observed  in  the  ray 
paths  (see  Figure  IV-C.5) .  These  changes  arise  mainly  from  the 
tidal  velocities  overcoming  the  group  velocities  during  part 
of  the  tidal  cycle.  This  causes  the  rays  to  propagate 
backwards  relative  to  their  group  velocity  resulting  in 
negative  Eulerian  frequencies;  only  the  rays  starting  close  to 
slack  water  at  low  tide  have  straight  paths.  This  behavior 
suggests  that  nonlinear  effects  become  dominant,  leading  to 
the  question  about  the  validity  of  the  plane  wave  assumption. 
Another  significant  difference  is  that  compared  with  the 
previous  runs,  the  reduced  vertical  group  velocities  do  not 
allow  the  initially  upward  propagating  rays  to  leave  the 
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upper,  strongly  stratified  layers,  causing  them  to  become 
trapped  near  the  surface.  Therefore,  independent  of  their 
validity,  the  upward  rays  are  neglected  in  the  rest  of  this 
discussion  because  of  their  negligible  effects  at  depths  where 
the  MBAY  observations  showed  strong  nonlinearities. 

The  reliability  of  the  ray  paths  for  run  #3  can  be 
evaluated  by  analyzing  the  behavior  of  variables,  such  as 
frequency  and  wavenumber,  and  by  evaluating  how  consistent 
they  are  with  the  known  tidal  flow  and  physics.  Figure  IV-C.6 
displays  the  time  series  of  wave  frequency  at  the  source.  The 
changes  in  the  Eulerian  frequency  are  dominated  by  doppler 
shifting  as  expected;  however,  from  the  intrinsic  frequency 
plot,  it  can  be  seen  that  the  straining  terms  also  introduce 
some  variability,  which  superimposed  with  the  doppler  shift 
due  to  the  advection  by  the  underlying  flow,  attenuates  the 
maximum  Eulerian  frequency  and  further  decreases  the 
frequency,  breaking  the  symmetry  of  the  doppler  shifting.  This 
effect  changes  the  mean  wave  frequency  at  the  source,  but  not 
significantly  enough  to  change  the  output  velocities.  However, 
this  effect  is  important  in  the  generation  of  residual 
currents,  together  with  the  amplitude  modulation  determined  by 
the  wave  action.  In  Figure  IV-C.7,  the  time  changes  of 
intrinsic  frequency  along  the  rays  can  be  seen.  Rays  starting 
during  the  flood  increase  their  intrinsic  frequency,  while 
rays  starting  during  the  ebb  show  a  decrease.  This  effect  is 
similar  to  those  discussed  for  runs  #1  and  #2  which  caused  ray 
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attenuation.  In  the  Euler ian  frequency  plot,  rays  starting 
during  the  flood  with  negative  frequencies  can  also  be  seen, 
indicating  that  they  are  traveling  in  the  opposite  direction 
to  their  group  velocity.  During  this  initial  period,  the  rays 
are  over  the  shelf  where  there  are  no  spatial  inhonogeneities 
in  the  velocity  field,  so  the  wavenumbers  will  not  change. 
This  can  be  seen  in  Figure  IV-C.8,  which  displays  the 
evolution  of  wavenumber  components  along  each  characteristic. 

To  compare  the  model  output  with  the  MBAY  data,  the 
velocity  field  at  each  instant  is  computed  using  the 
velocities  along  each  characteristic.  This  is  accomplished  by 
computing  the  intrinsic  frequency  at  each  point,  for  each 
time,  along  the  rays.  Based  on  these  values,  an  instantaneous 
velocity  amplitude  is  estimated  from  the  wave  action  value  at 
that  point.  The  individual  velocities  are  then  computed  using 
this  instantaneous  amplitude  and  intrinsic  frequency,  and  a 
phase  lag  determined  by  the  time  taken  to  arrive  at  that 
point,  as  shown  in  schematic  of  Figure  IV-C.9.  Also,  based  on 
the  interfacial  mode  for  a  step  slope  as  shown  in  Figure  IV- 
C.IO  (Baines,  1982),  the  initial  phase  at  the  shelf  break  of 
the  Internal  waves  propagating  offshore  is  estimated  to  be 
leading  the  tidal  forcing  by  w/2, 

A  47  by  100  velocity  grid  is  defined  enclosing  the 
continental  slope  between  the  shelf  break  and  9Km  offshore, 
and  from  the  surface  to  1000m  depth,  in  order  to  look  at  the 
resultant  velocity  field  of  a  ray  family.  Velocity  estimates 


along  each  ray  are  projected  on  to  the  x-z  plane  at  each 
instant.  The  velocity  field  between  two  consecutive  projected 
characteristics  is  computed  by  linear  interpolation  along  the 
z-dimension.  The  rays  are  computed  every  15min  through  a  tidal 
cycle  to  minimize  aliasing  effects,  and  to  limit  the  linear 
interpolation  needed  between  successive  rays.  However,  at 
greater  distances  from  the  shelf  break,  the  interpolation 
between  diverging  rays  produces  poor  velocity  estimates, 
resulting  in  broadening  of  some  of  the  features. 

To  study  the  sensitivity  of  the  resultant  output  flow 
field  as  a  function  of  the  number  of  rays,  the  density  of 
upward  rays  was  increased  to  start  every  7.5min  during  the 
tidal  cycle  for  run  #2.  No  significant  changes  were  found, 
indicating  that  the  ISmin  ray  spacing  is  appropriate  for  the 
grided  domain. 

A  24  hour  run  was  performed  to  estimate  residual  flows  for 
run  #3.  Figure  IV'-C.ll.a/b  show  maps  of  the  24  hour  averages 
of  U  and  V  velocity  components.  The  existence  of  the  velocity 
peaks  in  the  range  of  0.20m/s  to  O.lOm/s  in  the  velocity  field 
at  8250m  distance  between  30  and  70m  demonstrate  that 
nonlinear  interactions  of  the  internal  tides  with  the  tidal 
flow  are  able  to  produce  significant  residual  flows.  These 
depth  ranges  and  magnitudes  are  comparable  to  the  0.20m/s,  70m 
depth  baroclinic  northward  residual  velocities  observed  during 
the  MBAY  experiment  at  the  western  portion  of  the  surveyed 
area  (see  Chapter  III) .  These  residual  flows  are  associated 
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with  the  upper  portion  of  the  internal  wave  beam,  deepening 
further  away  from  the  shelf  break  and  showing  frequent 
reversals.  It  is  not  possible  to  verify  all  these  features 
from  the  NBAY  data,  where  the  northward  residual  flows  were 
observed  only  along  one  side  of  the  surveyed  box  preventing 
the  estimation  of  the  current  across-shore  extent.  However,  a 
reversal  in  direction  at  the  eastern  portion  was  observed  in 
the  form  of  a  narrow  southward  current  with  500- 1000m  width  as 
displayed  in  Figure  III-A.4. (ii) ,  centered  at  1000m  box  side 
distance.  This  flow  reversal  is  consistent  with  the  model 
results  that  predict  residual  flow  reversals  over  distances  of 
approximately  500m  in  the  cross-slope  direction. 

To  assess  the  way  in  which  these  residuals  occur,  time 
series  of  the  instantaneous  profile  velocity  components  were 
computed  at  three  different  locations  at  distances  650Gm, 
7500m,  and  8750m  along  the  slope  as  displayed  in  Figures  IV- 
12.a/b/c.  These  locations  were  chosen  to  be  representative  of 
far,  intermediate  and  close  fields  from  the  shelf  break  wave 
source.  The  profiles  at  6500m  shows  wavelike  oscillations  with 
the  U  and  V  components  close  to  quadrature.  These  oscillations 
change  frequency  with  depth,  such  that  a  24  hour  periodicity 
can  be  observed  in  the  V  component  at  70m  depth,  although  the 
U  component  at  the  same  depth  has  a  12  hour  periodicity.  At 
greater  depth,  higher  frequencies  are  observed  with  a  6  hour 
period  defined  in  the  layers  enclosed  by  the  wave  beam.  This 
depth  dependent  frequency  is  consistent  with  the  MBAY  data. 
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which  shows  24  hour  periods  in  the  V  component  in  areas  and 
depth  ranges  where  current  residuals  occur;  however,  there  is 
no  evidence  for  these  higher  frequency  waves  in  the  data.  This 
may  be  explained  by  the  two  hour  sampling  strategy  used  that 
might  have  aliased  them,  or  smearing  of  these  features  by 
effects  such  as  a  broadband  baroclinity  tidal  source  and  real 
variations  in  topography.  The  profiles  at  7500  and  8750  do  not 
show  the  low  frequency  component;  however,  they  display  a 
strong  frequency  modulation  at  the  deeper  portion  of  the  beam 
with  periods  decreasing  from  semidiurnal  to  a  few  hours. 
Therefore,  as  in  the  NBAY  observations,  the  diurnal 
oscillations  of  the  modeled  velocity  field  are  spatially 
limited.  Comparing  the  profiles  at  7500m  and  8750m,  it  can  be 
observed  that  the  higher  frequencies  become  stronger  near  the 
source,  indicating  that  this  area  will  be  more  sensitive  to 
model  parameters  values  and  to  instabilities. 

The  profile  time  series  show  that  in  the  far  field, 
currents  within  the  beam  are  nearly  constant  with  depth,  as 
Figure  IV-C.12.a  shows;  however,  strong  vertical  shear  is 
observed  in  the  upper  layers,  close  to  depths  where  the  lower 
frequency  oscillations  are  observed.  On  the  other  hand,  the 
intermediate  and  near  field  velocities  displayed  in  Figures 
IV~C.12.b/c  become  strongly  sheared  after  24  hours  as  strong 
vertical  velocity  changes  are  observed  in  the  upper  layers.  To 
assess  potential  shear  Instabilities  that  might  force  mixing 
within  the  internal  wave  beam,  a  Richardson  number  criteria 
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was  calculated  using  the  hourly  velocity  fields  over  the  slope 
and  the  stratification  displayed  in  Table  IV-C.l.  The  number 
of  occurrences  of  critical  or  below  critical  conditions  of 
Rj<0,25  were  counted  as  displayed  in  Figure  IV-C.13.  The  high 
occurrence  of  critical  conditions  shows  the  potential  of  this 
mechanism  to  force  shear  instabilities  similar  to  those 
observed  during  the  MBAY  experiment;  however,  this  result 
should  change  significantly  with  geometry  and  barotropic  tidal 
flow  choices,  so  no  direct  correspondence  would  be  expected 
between  observations  and  this  model  result. 

These  model  results,  which  reproduce  features  of  the  MBAY 
observations,  show  that  nonlinear  interactions  of  the  internal 
wave  field  generated  at  the  shelf  break  with  the  barotropic 
tide  are  candidate  mechanisms  for  the  rectification  of 
internal  tides.  However,  the  modeled  velocity  field  does  not 
totally  explain  the  spatial  structure  of  the  observed  data, 
though  they  occur  at  similar  depth  ranges  and  with  similar 
magnitudes.  This  could  be  due  to  the  existence  of  enhanced 
internal  waves  activity  associated  with  the  canyon,  and  due  to 
greatly  simplified  model  geometry  and  barotropic  tidal  flow. 
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V.  CONCLU8IOMS 


A  process  study  over  a  coastal  shelf  break  was  developed 
to  examine  the  role  of  tidal  dynamics  in  mixing  under 
conditions  of  steep  topography  and  strong  stratification.  A  3 
day  survey  (MBAY)  of  the  upper  ocean  was  conducted  close  to 
the  mouth  of  Monterey  canyon  measuring  velocity,  temperature 
and  salinity  profiles  along  a  4km  square  box.  Pycnocline 
displacements  of  20m  with  wavelengths  above  10000m  were 
observed  and  discussed  in  Chapter  III-B.  Chapter  IV-A  showed 
that  they  result  from  an  internal  tide  generation  process, 
which  also  forced  internal  waves  below  the  pycnocline.  Strong, 
rectified  bursts  of  shear  kinetic  energy  were  observed  at 
depths  from  50  to  100m  with  approximately  24  hour  periods.  The 
discussion  in  Chapter  III-B  showed  that  the  measured  weak 
winds  cannot  force  either  the  observed  highly  inhomogeneous 
and  strong  current  field  and  vertical  shear,  or  the  observed 
semidiurnal  thermocline  oscillations.  Using  the  studies  of  the 
effects  of  geostrophic  flow  over  a  canyon  by  Klinck  (1988, 
1989)  and  the  data  of  Koehler  (1990),  the  cross-canyon 
geostrophic  flow  is  expected  to  be  much  less  energetic  than 
the  tidal  flow.  This  is  confirmed  by  the  previous  work  in  the 
area  siunmarized  in  Breaker  et  al.  (1989)  and  by  the 
unpublished  current  meter  buoy  data  kindly  made  available  by 
C.H.Pilskaln,  (Pillsbury,  et  al.,  1992).  Therefore,  the 
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effects  of  geostrophic  ilow  across  the  canyon  on  the 
generation  of  waves  and  residual  flows  cannot  account  for 
these  observations.  It  was  shown  in  Chapter  IV-B  that  the 
observed  bursts  of  energy  could  result  from  the  rectification 
of  the  tidal ly  forced  internal  waves  undergoing  interaction 
with  the  barotropic  tidal  flow.  An  Ir.^'^rnal  mixed  layer  at  the 
thermocline  depth  with  40m  depth  is  also  observed  and 
discussed  in  Chapter  III-C.  It  was  founu  that  this  mixed  layer 
was  probably  forced  by  shear  supplied  by  the  non-linear 
internal  waves  below  the  thermocline,  and  by  shear  associated 
with  the  thermocline  oscillations.  However,  the  data  and  the 
theoretical  analysis  did  not  allow  a  positive  identification 
of  the  triggering  mechanism. 

The  MBAY  data  sets  and  the  numerical  results  presented  in 
chapter  IV  show  the  need  for  a  complete  analysis  of  the 
tidal ly  forced  dynamics  of  the  Monterey  Bay  as  suggested  by 
Shea  and  Broenkow  (1982).  The  good  comparison  of  observed 
thermocline  oscillation  amplitude  and  frequency  with  the 
linear  internal  tide  generation  model  by  Baines  (1982)  showed 
continental  slope  tidal  dynamics  to  be  valid  near  the  canyon 
mouth.  Furthermore,  velocity  observations  taken  at  the  canyon 
mouth  were  compared  with  the  numerical  velocity  fields 
resulting  from  a  nonlinear  internal  wave  propagation  model, 
using  model  parameters  typical  of  those  in  Monterey  Bay,  but 
considering  a  linear  continental  slope.  The  model  results 
suggest  the  existence  of  multiple  wave  frequencies  and 


92 


nonlinearities  resulting  in  strong,  unidirectional,  spatially 
inhonogeneous  currents,  as  observed  during  the  MBAY 
experiaent.  However,  there  is  no  direct  correspondence  between 
observations  and  model  results  because  of  the  simplified 
geometry  and  barotropic  tidal  model.  It  was  also  shown  in  the 
present  work  that  internal  tide  dynamics  in  the  area  are 
largely  dominated  by  nonlinear  phenomena  that  might  be 
enhanced  by  the  actual  steeper,  more  irregular  topography  and 
inhomogeneous  stratification  which  are  found  inside  Monterey 
canyon.  The  current  approach  can  be  extended  to  incorporate 
these  cases;  however,  better  internal  wave  dissipation  models 
and  higher  order  corrections  should  be  developed  and  studied 
to  increase  the  confidence  of  the  numerical  results.  A  more 
detailed  discussion  of  these  conclusions  follows. 

A.  ZMTKtIIAL  TZDB8  AT  THE  MODTB  OF  THE  NOMTEREY  CAMYOM 

Based  on  the  work  by  Baines,  (1983),  it  can  be  assumed 
that  internal  tide  dynamics  at  the  canyon  mouth  are  similar  to 
those  over  the  open  shelf,  as  discussed  in  Chapter 
The  effects  of  geostrophic  flow  over  canyons  has  been  studied 
by  Klinck  (1988; 1989),  and  based  on  these  papers  and  the 
current  mooring  data,  the  effects  of  geostrophic  flow  (i.e., 
eddies  from  the  California  Current)  over  the  canyon  can  be 
neglected  for  the  MBAY  case. 

It  was  shown  in  Chapter  IV-A  that  the  estimated  barotropic 
tidal  velocity  field  for  a  strong  thermocline,  as  observed 
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during  MBAY,  will  produce  a  maximum  internal  tide  forcing 
close  to  the  thermocline  at  the  shelf  break.  Furthermore, 
using  the  generation  model  by  Baines  (1982) ,  it  was  found  in 
Chapter  IV**A  that  there  will  be  an  interfacial  mode  of  O(lOm) 
amplitude  and  0(20Km)  wavelength  at  the  thermocline.  This  mode 
decouples  the  upper  layer  from  the  internal  wave  propagation 
in  the  stratified  continuum  below,  and  forces  shear  at 
thermocline  depth  that  produces  Richardson  numbers  close  to  1 
during  parts  of  the  tidal  cycle.  These  results  are  in 
agreement  with  the  measurements  discussed  in  III-B,  before 
hour  152  in  the  West  portion  of  the  surveyed  area  where 
semidiurnal  thermocline  oscillations  with  10m  amplitudes  and 
wavelengths  above  lOKm  where  observed.  Along  the  east  leg, 
close  to  the  canyon  southern  wall,  the  observations  did  not 
show  the  same  tidal  dependence.  This  suggests  that  other 
effects  at  the  canyon  wall  are  acting,  possibly  by  canceling 
internal  tide  mechanisms  or  by  creating  a  shadow  £one.  These 
mechanisms  might  be  related  to  barotropic  velocities  across 
the  canyon  (alongslope) .  Unfortunately,  the  observations  and 
the  theoretical  formulation  do  not  offer  any  conclusions  about 
this  case. 

B.  BAROCLIMIC  RECTIFICATION  OF  TIDES 

Barotropic  tidal  rectifications  over  the  continental  slope 
close  to  Monterey  can  occur  by  the  generation  of  residual 
vorticity  due  to  the  conservation  of  potential  vorticity 
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(dominated  by  vortex  stretching) ,  as  discussed  in  Chapter  III- 
D.2.  Estimates  for  residual  vorticity  (Robinson,  1982)  in  the 
present  case  are  of  the  order  3.6xl0'^s'\  resulting  in 
northward  alongslope  barotropic  velocities  of  about  0.05ms'\ 
close  to  the  shelf  break.  These  estimates  compare  well  with 
the  observations  as  mentioned  in  Chapter  III-D.2.  However,  the 
NBAY  data  shows  that  the  vertical  mean  values  are  a  result  of 
vertically  layered  residuals  of  jet-like  current  pulses  which 
occur  at  24  hours  intervals.  The  non-linear  modal  approach  for 
a  stratified  ocean  (following  Mass  and  Zimmermann,  1989a, b) 
showed  that,  although  there  is  much  in  common  with  the 
homogeneous  case,  the  main  effect  of  including  vertical 
stratification  in  tidal  rectification  by  vortex  stretching  is 
a  bottom  intensification,  which  cannot  be  detected  in  the 
near-surface  MBAY  observations  and  is  not  consistent  with  the 
current  maxima  observed  within  the  water  column.  Nonlinear 
layer  models  for  the  MBAY  case  can  predict  residual  currents 
in  the  upper  layer  of  about  0.02ms'^  and  0.002ms''  in  the  lower 
layer,  values  that  are  too  small  compared  with  the 
observations.  However,  based  on  Baines  (1982),  the  velocity 
amplitudes  of  the  internal  waves  forced  over  the  shelf  break 
were  estimated  to  be  of  about  0.15ms'',  resulting  in  estimated 
residuals  up  to  0.05ms''.  Again,  this  value  is  smaller  than  the 
observed  layered  time  residuals  which  are  approximately  0.25ms' 
'.  To  evaluate  the  nonlinear  processes  that  could  lead  to 
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these  rectified  velocities,  a  ray  tracing  model  was  developed 
by  considering  the  effects  of  advection  and  straining  of  the 
internal  wave  field  by  the  barotropic  topographically  enhanced 
tidal  flow.  The  results  of  this  model  showed,  for  a  simplified 
continental  slope  similar  to  that  in  Monterey  Bay,  24  hour 
velocity  residuals  between  0.20ms'^  to  0.10ms*\  at  40m  to  70m 
depths  along  the  upper  portion  of  an  internal  wave  beam 
starting  at  thermocline  depth  (~25m)  at  the  shelf  break.  These 
rectified  velocities  start  at  distances  as  close  as  IKm  from 
the  shelf  break,  and  can  be  observed  as  far  as  4Km  from  the 
source,  showing  sign  reversals  within  distances  as  short  as 
500m.  These  residuals  have  the  same  order  of  magnitude  as  the 
KBAY  observed  values  of  0.25ms'\  and  the  depth  ranges  are 
comparable.  However,  for  the  MBAY  data  a  determination  of  the 
acrossoshore  extent  of  the  observed  strong  northward  velocity 
residues  is  not  possible,  though  Figure  III-A.4.(ii)  shows  a 
well  defined  residual  velocity  of  0.15ms'^  in  the  opposite 
direction,  extending  over  a  cross~shore  distance  of  500**1000m 
in  agreement  with  the  model  results.  However,  the  depth  ranges 
over  which  the  residuals  are  observed  are  broader  than  the 
model  results.  This  is  probably  due  to  many  effects,  such  as 
the  existence  of  pressure  gradients  established  in  the  canyon, 
which  can  enhance  the  internal  wave  field  (Baines,  1983; 
Huthnance,  1989) ,  and  due  to  the  simplified  tidal  flow  and 
topography  which  were  used  in  the  model. 
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The  tenporal  variability  of  the  model  currents  along  the 
wave  bean  show  a  nultiplicity  of  frequencies  with  smaller 
values  close  to  diurnal  frequencies  and  at  locations  where  the 
maximum  current  residuals  were  observed;  on  the  other  hand, 
higher  frequencies  were  observed  close  to  the  source  at  the 
shelf  break.  This  result  partially  agrees  with  the  MBAY 
observations,  which  showed  residual  velocities  resulting  from 
diurnal  jet-like  features  at  depths  between  50  to  100m,  while 
semidiurnal  frequencies  were  observed  at  other  depths. 
However,  the  MBAY  observations  showed  the  low  frequency 
fluctuations  as  deep  as  I50m  throughout  the  layers  where  the 
residuals  were  observed.  Again,  this  discrepancy  between  the 
model  and  observations  is  likely  due  to  the  existence  of  the 
canyon  and  due  to  the  simplified  topography  and  stratification 
used  in  the  model.  The  MBAY  data  did  not  resolve  higher 
frequency  oscillations  due  to  the  sampling  strategy.  Therefore 
it  is  not  possible  to  make  a  complete  comparison  with  the 
model  results. 

Based  on  this  discussion,  nonlinear  interactions  of  the 
internal  waves  generated  at  the  shelf  break,  together  with  the 
barotropic  and  topographically  changed  tidal  flow,  are  a 
mechanism  that  can  cause  the  occurrence  of  diurnal,  strong, 
baroclinic,  unidirectional  jet  like  features,  similar  to  those 
observed  during  the  MBAY  experiment.  These  nonlinear 
interactions  can  also  contribute  to  current  shear  residuals 
with  the  current  reversing  within  small  vertical  distances. 
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C.  UrTBSMAL  NXZXMG  DUB  TO  THE  TXDALLY  FORCED  SHEAR 

The  nonlinear  internal  tide  model  developed  in  Chapter  IV, 
showed  multiple  frequencies  and  a  strongly  baroclinic  current 
structure,  occurring  along  a  beam  which  starts  at  the  shelf 
break.  The  modeled  vertical  current  shear  was  obsezrved  to  be 
stronger  near  the  source,  where  higher  frequencies  exist,  and 
at  each  location  close  to  the  upper  limit  of  the  wave  beam, 
where  the  lower  frequencies  were  observed.  A  Richardson  number 
criteria  was  applied  to  the  modeled  internal  tide  field  to 
estimate  regions  where  mixing  was  likely  to  occur.  The  number 
of  occurrences  of  Richardson  number  below  the  critical  value 
of  0.25  were  counted  for  a  24  hour  model  run.  It  was  found 
that  internal  tide,  wave-induced  shear  instability  mechanism 
could  account  for  the  occurrence  of  local  internal  shear 
instabilities,  which  will  dissipate  wave  energy  and  force 
internal  mixing  layers,  similar  to  the  observation  during  the 
MBAY  experiment  within  the  thermocline. 

D.  FUTURE  WORK 

Assumptions  in  the  numerical  model  developed  in  Chapter  IV 
are  not  completely  satisfied  for  the  case  under  consideration. 
The  comparisons  with  the  MBAY  observations  do  not  completely 
validate  the  numerical  results.  Therefore,  it  is  desirable  to 
complement  this  work  by  estimating  the  magnitudes  of  higher 
order  corrections,  and  by  comparing  results  with  data  other 
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than  the  MBAY  survey  or  with  other  existing  validated  model 
results. 

When  comparing  the  numerical  results  with  observations, 
significant  discrepancies  were  found  and  attributed  to 
simplified  geometry  and  barotropic  tidal  flow.  Therefore,  one 
obvious  extension  of  the  present  work  is  to  improve  the 
barotropic  tide  model  and  to  use  real  topography,  including 
the  effect  of  the  canyon,  using  an  approach  similar  to  the  one 
developed  by  Baines  (1983). 

Other  improvements  of  the  numerical  model  can  be  achieved 
by  including  wave-wave  interactions,  better  identification  and 
handling  of  some  wave-current  interactions  where  the  ray 
tracing  method  is  clearly  not  valid,  and  by  using  a  more 
realistic  model  of  the  surface  and  bottom  reflections 
(currently  they  are  assumed  to  be  perfectly  elastic) . 
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APPBMDZX  A  -  TABLES 
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Tabl«  Xl>A.l 


Geographycal  coordinates  of  the  surveyed  square  area  during 
the  MBAY  experiment. 


LATITUDE 

LONGITUDE 

Corner  #  1  (NE) 

36.7366 

122.0033 

Comer  I  2  (SE) 

36.7022 

122.0033 

Corner  #  3  (SW) 

36.7022 

122.0483 

Corner  #  4  (NW) 

36.7366 

122.0483 

Center  (P.Lobos) 

36.72 

122.025 
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TABLE  IV-A.l 


Suanuiry  of  significant  parameters  of  Baines  (1982)  model, 
showing  numerical  values  estimated  for  MBAY  conditions. 
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DBFZlIXTZON 


2.ir/Tidal 


MBAY  VALUES 


1 . 4*10‘^S'’ 


2.9xl0‘^s‘’-  Mean 


4xl0‘®s‘’-  Sfc 


20xl0'®s'’-  2  On 


3xl0'^s*’-  lOOn 


2x10V'-  lOOOn 


0.03  -  Sfc 


Tidal 


Frequency 


Brunt-Vaisala 

Frequency 


Internal  Wave 


Slope 


Paraneters 


determining 
the  existence 
of  an 

interfacial 


mode 
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TMU  ZT-O.l 


SuMBary  of  model  runs  used  to  dleouee  the  nuaerloal  procedure 
and  to  eetlaate  the  flow  aaeooiated  with  nonlinear  effecta  on 
Internal  waves  propagating  frcm  the  shelf  break. 


RUNfl 

RUN«2 

RUN  #3 

Slope 

Attenuatted 

MBAY 

MBAY 

Hor.  Scales 

**2  Km 

~7Km 

~7  Km 

Period 

~2Hr 

-6Hr 

12  Hr 

Stratification 

MBAY 

MBAY 

MBAY 

Note:  MBAY  Stratification  ccHTespoads  to  the  foliowing  values  of  N: 


H 

0 

4X10** 

20 

20x10** 

30 

2.9X10** 

100 

2.9X10** 

1000 

2X10** 
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APPEMDZX  B  -  FIGURES 


Figur*  I-A.l 


Physical  processes  assumed  Important  in  a  coastal  region  with 
a  highly  irregular  topography.  The  processes  inside  the 
thicker  line  apply  to  the  Monterey  Bay  Mixing  Experiment 
(MBAY) ,  in  September  1991. 
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rigur*  zz-A.i 


Figure  (a)  shows  a  nap  of  the  Monterey  Bay.  The  Monterey  Bay 
Mixing  Experiment  (MBAY)  location  and  the  mooring  site  of  the 
buoy  data  of  C.H.Pilskaln  (Pillsbury  et  al.,  1992)  are  marked 
with  arrows.  Figure  (b)  shows  a  schematic  of  the  MBAY 
experiment  and  a  zoom  of  ocean  floor  topography.  The  RV  "Pt. 
Sur"  performed  repeated  profiling  around  a  square  path 
approximately  every  2  hours  during  three  days  acquiring  ADCP 
velocity  profiles,  586  Tow^Yo  CTD  profiles  and  meteorological 
data.  Concurrently,  the  RV  "Pt.  Lobos"  was  stationed  close  to 
the  center  of  the  box  controlling  a  ROV  performing  direct 
turbulence  measurements.  The  surveyed  area  covered  a  4  Xm 
square  in  the  South-West  wall  of  the  Monterey  canyon,  whose 
coordinates  are  shown  in  Table  IIl-A-1.  This  area  has  a  highly 
Irregular  topography  with  depth  ranging  from  700m  at  the 
north-west  comer  to  1500m  at  the  north  and  south  legs. 
Velocity  data  were  gathered  down  to  250m  depth  and  ancillary 
data  (temperature  and  salinity)  were  observed  to  80m  depth 
along  the  surveyed  box.  Strong  shear  and  an  internal  mixed 
layer  were  observed  at  the  western  side  of  the  box. 
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Mooring  site 


Flgurn  II-A.l  (a) 


MBAY 

MIMNOEXPERIIIENT 


MonlartyBay 


North  Log 


TOPOGRAPHY 
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rigor*  ii-A.l  (b) 


rlgur*  11-A.2 


Along  track  time  series  of  eastward  (U)  and  northward  (V)  ADCP 
velocities  from  7  to  120  meter  depth  for  box  #25.  Each 
vertical  profile  has  been  demeaned  relative  to  a  reference 
layer  at  100  meter.  Profiles  were  gathered  every  1  minute.  The 
profiles  corresponding  to  each  box  corner  are  marked  along  the 
x-axls.  Strong  baroclinic,  jet  like  features  are  evident  in 
both  components.  These  jets  appear  to  be  decouple  from  the 
surface  and  to  reach  maximum  values  near  the  thermocline.  The 
strong  shear  layers  at  approximately  20m  and  50m  depths 
between  ADCP  profiles  4050  and  4070  (better  seen  in  the  U 
component  profiles) ,  and  are  simultaneous  with  the  isotherms 
displacements  (upward  at  20m  and  downward  at  50m)  shown  in  the 
along  track  temperature  profiles  (between  Tow-Yo  profiles  472 
and  481)  in  Figure  III-A.3. 
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Flgur*  II-A.3 


Along  track  time  series  of  (a)  sigma-theta  and  temperature, 
and  (b)  salinity,  from  0  to  80  meter  depth  for  box  #25.  The 
profiles  at  each  box  corner  are  marked  along  the  x-axis.  The 
deepening  of  isotherms  appears  to  be  related  to  the  occurrence 
of  the  strong  shear  shown  in  Figure  III-A.2.  There  is  some 
evidence  of  temperature  inversions  that  are  being  compensated 
by  salinity.  Internal  waves  are  evident  at  the  thermocline  in 
the  sigma-theta  and  temperature  profiles. 
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Figure  ZZ-A.3  <b) 


116 


Figur*  ZI-A.4 


T/S  diagrams  of  the  along  track  Tow-Yo  data  over  the  West  most 
leg  12  hours  apart:  (a)  corresponds  to  box  #1  at  01:00  in 
9/5/91;  and  (b)  corresponds  to  box  #5  at  12:00  of  the  same 
day.  The  new  water  type,  appearing  at  the  depth  of  30  meter  in 
(b)  (box  #5) ,  corresponds  to  a  loss  in  potential  energy 
(denser  water)  and  must  be  a  result  of  mixing  or  intrusions. 
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PLOT  FOR  MBRYS  -  BOX  #1 


(nsd) 


8  PLOT  FOR  MBRYS  *BOX#S 


rigur«  iz*-A.4  (b) 
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(nsd 


rigur* 


Figure  (i)  shows  the  2-dimensional  autocorrelation  function 
used  for  the  assimilation  of  the  MBAY  data  into  a  grid.  Figure 
(ii)  shows  the  correlation  values  at  a  zero  distance  lag  as  a 
function  of  the  time  lag.  Similarly,  Figure  (iii)  shows  the 
correlation  values  at  a  zero  time  lag  as  a  function  of  the 
special  lag.  Correlation  length  and  time  scales  were  chosen 
based  on  the  sampling  scale.  The  scaled  variables  'tt'  and 
'rr*  correspond  to  the  squares  of  the  scaled  time  and  distance 
lags  and  are  displayed  in  Figures  (iv)  and  (v) . 
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Tignxm 
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Figur*  ll-B-1.2 


An  example  of  the  implementation  of  the  optimum  interpolation 
scheme  used  to  analyze  the  several  observations  onto  a  grid. 
An  empirical  correlation  function  with  450  meter  and  15 
minutes  scales  was  used.  The  upper  Figure  (i)  shows  the  depth 
averaged  number  of  data  points  used  for  each  of  the  27  runs, 
around  the  surveyed  box.  Less  than  5  data  points  were  used  for 
most  of  the  grid  points  along  the  boxes.  The  lower  Figure  (ii) 
shows  contours  of  the  depth  averaged,  normalized  KMS  errors, 
at  each  grid  point  for  each  box.  RMS  errors  are  indicators  of 
the  dominant  distances  between  the  data  points  and  each  grid 
point.  Gaps  in  the  data  were  assigned  RMS  of  1. 
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rigur«  ll-B-1.2 
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rlgur* 


Absolute  U  (i)  and  V  (ii)  velocity  conponent  spectra 
(unfiltered  data) ,  after  applying  the  GPS  navigated  ships 
velocity.  Three  bands  are  identified.  The  higher  frequency 
band  is  assumed  to  be  instrumental  noise  and  was  filtered  from 
the  data  using  a  4th  order,  zero-phase,  Butterworth  filter, 
within  the  error  band,  the  observed  variance  is  assximed  to  be 
from  instrumental  noise  and  from  physical  current  stinicture. 
This  band  was  used  to  estimate  the  error  of  the  absolute 
current  estimates  (app.  4  cm/s  at  95%  confidence) .  The 
variance  in  the  structure  band  is  assumed  to  dominated  by 
physical  current  structure.  The  peak  at  approximately  2  hour 
period  results  from  the  strong  inhomogeneities  in  the  mean 
velocity  field  over  the  sampled  area.  The  semi-diurnal  tidal 
component  is  observed  in  both  U  and  V  spectra.  The  peaks  at 
approximately  4  hour  period  have  no  obvious  explanation.  They 
might  be  due  to  doppler  shifts  of  the  tides  or  to 
inhomogeneous  spatial  current  structure. 
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Figure  ZII-A.l 


Vector  plot  of  the  time  series  of  winds  (i),  during  the  MBAY 
experiment.  The  wind  magnitudes  were  low  to  moderate  with  a 
rapid  change  in  direction  at  hour  156.  The  magnitude  and 
direction  of  the  wind  are  displayed  in  (ii) .  Below  the  figure, 
the  hours  of  the  predicted  high  and  low  tides  are  marked  along 
the  x-axis. 
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Tigura  izz-A.l 
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Figur*  III-A.2 


Vector  plot  of  vertically  averaged  (upper  250ia)  absolute 
velocities  averaged  over  all  27  box  measurements.  Stronger 
currents  are  observed  at  the  west  leg  over  the  Southward 
facing  slope  of  the  canyon  (Figure  II-A.l).  Overall,  an 
anticyclonic  circulation  can  be  seen.  Figure  III-A.6  gives 
quantitative  estimates  of  this  circulation  and  of  the 
divergence . 
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MBAY  -  TIME  AND  VERTICALLY  AVERAGED  VELOCITY  FIELD 


Figur*  ZXI-JL.2 
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East- West  axis  Distance  In  meters 


rigur*  IZZ-A.3 


Profiles  of  time  averaged  eastward  (U)  velocity  component  for 
the  3  day  measurements.  Each  panel  represents  the  observed 
values  along  one  leg.  Figure  (i)  corresponds  to  the  east  most 
leg,  (ii)  corresponds  to  the  south  most  leg,  (iii)  corresponds 
to  the  west  most  leg  and  finally  (iv)  represents  the  north 
most  leg.  Each  panel  starts  at  the  beginning  of  the 
corresponding  leg  in  the  following  direction  (NE-SE-SW-NW-NE) . 
They  show  a  strong  baroclinic  and  spatial  inhomogeneous  mean 
flow.  See  the  text  for  other  figure  interpretations. 
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Figure  III-A.3 


131 


Figur«  XI1-A.4 


Profiles  of  time  averaged,  northward  (V)  velocity  component 
for  the  3  day  measurements.  Each  panel  represents  the  observed 
values  along  one  leg.  Figure  (i)  corresponds  to  the  east  most 
leg,  (ii)  corresponds  to  the  south  most  leg,  (iii)  corresponds 
to  the  west  most  leg  and  (iv)  represents  the  north  most  leg. 
Each  panel  starts  at  the  beginning  of  the  correspondent  leg, 
in  the  following  direction  (NE-SE-SW-NW-NE) .  They  show  a 
strong,  baroclinic  and  spatial  inhomogeneous  mean  flow.  See 
the  text  for  other  figure  interpretations. 
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Figur*  III-X.5 


Time  series  of  box  mean  divergence  times  the  inertial 
frequency  (a) ,  box  mean  vorticity  (b)  and  time  derivative  of 
the  box  mean  vorticity  (c) .  The  divergence  was  computed  based 
on  the  normal  component  of  the  current  relative  to  the  ship's 
closed  track  and  the  vorticity  based  on  the  along>track 
velocity  component.  The  tidal  signature  can  be  seen  on  the 
vertical  averages  of  the  divergence  (suggesting  either 
internal  wave  motion  or  topographically  forced  vertical 
velocities) .  A  strong  baroclinicity  is  evident  in  all  current 
components.  One  should  also  note  the  strong,  rapid  vorticity 
change  in  the  vertical  at  30m  depth  at  hour  152.  Section  III-C 
discusses  likely  consequences  of  this  vorticity  pattern.  The 
time  derivative  of  vorticity,  though  it  shows  the  correct 
order  of  magnitude  relative  to  the  divergence  that  one  can 
expect  from  linear  theory  for  an  homogeneous  ocean,  does  not 
follow  closely  the  divergence  patterns.  This  is  likely  due  to 
the  presence  of  near-inertial  internal  waves.  The  tidal 
signature  is  also  apparent  in  the  vertical  averages  but  is 
weaker  than  the  divergence  counterpart. 
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Box  Mean  Divergence  Estimates  *  corioiis 


rigor*  ZZZ-A.5  (a) 
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Figur«  IZZ-A.5  (b) 
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Figure  IIZ'A.5  (o) 
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Figux* 


Profiles  of  the  tenporal  evolution  of  the  absolute  value  of 
the  current  speed  at  4  points  selected  along  the  west  most 
leg.  The  hours  shown  in  the  x^axis  are  referred  to  the  0  hour 
of  September  1  1991.  Baroclinic  bursts  of  absolute  velocity 
are  observed  with  an  apparent  period  of  24  hours.  See  the  text 
for  a  more  detailed  interpretation  of  the  figure. 
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Flgur*  IZZ-B.! 


Contours  of  the  vertical  structure  of  tenperature  at  hour  124 . 
Figure  (i)  corresponds  to  the  east  most  leg,  (ii)  corresponds 
to  the  south  most  leg,  (iii)  corresponds  to  the  west  most  leg 
and  (iv)  represents  the  north  most  leg.  Temperature  data 
between  the  surface  and  8m  depth  corresponds  to  a  linear 
interpolation  between  the  sea  surface  temperature  and  the 
first  CTD  Tow-Yo  data  bin.  Each  panel  represents  the  observed 
values  along  one  leg.  Each  panel  starts  at  the  beginning  of 
the  correspondent  leg  in  the  following  directions  (NE-SE-SH- 
NW'HE) .  See  the  text  for  further  interpretations. 


East  Most  Leg  at  124  Hour  (  i  )  South  Most  Leg  at  1 24  Hour  (  |7  ) 


Figure  ZII-B.l 
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rlgur* 


Contours  of  the  vertical  structure  of  temperature  at  hour  130. 
Figure  (i)  corresponds  to  the  east  most  leg,  (ii)  corresponds 
to  the  south  most  leg,  (iii)  corresponds  to  the  west  most  leg 
and  (iv)  represents  the  north  most  leg.  Temperature  data 
between  the  surface  and  8m  depth  corresponds  to  a  linear 
interpolation  between  \.Iie  sea  surface  temperature  and  the 
first  CTD  Tow“Yo  data  bin.  Each  panel  represents  the  observed 
values  along  one  leg.  Each  panel  starts  at  the  beginning  of 
the  correspondent  leg  in  the  following  directions  (NE-SE-SW- 
NW-NE) .  See  the  text  for  further  interpretations. 
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Figure  III-B.3 


Velocity  contours  showing  the  time  evolution  around  the 
seunpled  box  at  19m  depth  ( i .  e . ,  approximately  the  mean 
thermocline) .  Figures  (a)  and  (b)  correspond  to  the  U  and  V 
current  components  respectively.  Figure  (b)  shows  northward 
currents  occurring  along  the  west  leg  with  a  tidal 
periodicity,  though  at  hours  125,  170  and  185  some  current 
enhancement  can  be  observed.  These  currents  are  not  similar  to 
those  shown  in  Figure  III-B.4.b,  that  show  a  24  hour 
periodicity;  however,  the  northward  strong  surface  currents 
between  hours  180-185  in  Figure  IV-B.3.b  are  consistent  with 
these  observations.  One  should  note  that  between  hours  140  and 
155  the  weaker  currents  at  19m  are  almost  simultaneous  with 
the  occurrence  of  strong  vertical  shear  at  this  depth  along 
the  west  leg,  as  can  be  seen  in  Figure  III-C.3,  suggesting 
that  part  of  the  kinetic  energy  might  have  been  used  for 
mixing  due  to  a  shear  instability  process. 
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Figur*  III-B.3  (b) 
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Figur*  III-B.4 


Velocity  component  contours  showing  the  time  evolution  around 
the  sampled  box  at  11m  depth.  Figures  (a)  and  (b)  correspond 
to  the  U  and  V  current  components  respectively.  Figure  (a) 
shows  a  near  surface  eastward  current  component  throughout  the 
box  until  hour  156,  when  this  component  shifts  to  a  westward 
direction  throughout  the  box.  This  shift  seems  to  propagate 
along  the  west  leg  between  hours  150  and  160.  Figure  (b)  show 
that  the  V  component  at  this  depth  is  following  the  tidal 
cycles  until  hour  165,  after  which  a  strong  surface  northward 
current  starts  at  the  beginning  of  the  west  leg  and  spreads 
toward  the  north-west  corner  between  hours  165  and  180.  After 
hour  180  the  current  reaches  a  maximum  and  becomes  spatially 
confined  to  the  western  part  of  the  surveyed  area. 
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ABSOLUTE  CURRENT  -  U  Component 
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Figure  XII-B.4  (b) 
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Figttr*  ZZI-B.5 


Velocity  contours  showing  the  tine  evolution  around  the 
sampled  box  at  51m  depth.  Figures  (a)  and  (b)  correspond  to 
the  U  and  V  current  components  respectively.  Figure  (a)  shows 
a  U  component  following  the  tidal  cycle  until  hour  170.  Rapid 
temporal  changes  occur  at  the  north-west  corner  after  hour 
150.  Figure  (b)  shows  the  a  strong  northward  V  component  along 
the  west  leg  corresponding  to  bursts  of  kinetic  energy  24 
hours  apart,  in  agreement  with  Figure  IV-A.5.  At  the  end  of 
the  time  series,  the  following  burst  is  weaker  than  the 
previous  two.  This  attenuation  is  simultaneous  with  the 
occurrence  of  stronger  currents  at  the  surface  of  Figure  IV- 
B.3.b. 
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Figur*  ZXI-B.5  (») 
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rigur*  XII-B.5  (b) 
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Flgur*  III-C.l 


Contours  of  the  vertical  structure  of  temperature  at  hour  168. 
Figure  (i)  corresponds  to  the  east  leg,  (ii)  corresponds  to 
the  south  leg,  (iii)  corresponds  to  the  west  leg  and  (iv) 
represents  the  north  leg.  Temperature  data  between  the  surface 
and  8m  depth  is  a  linear  interpolation  between  the  sea  surface 
temperature  and  the  first  CTD  Tow-Yo  data  bin.  Each  panel 
represents  the  observed  values  along  one  leg.  Each  panel 
starts  at  the  beginning  of  the  correspondent  leg  in  the 
following  directions  (NE-SE-SW-NW-NE) .  See  the  text  for 
further  interpretations. 
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Figure  III-C.2 


Contours  of  the  vertical  structure  of  temperature  at  hour  174. 
Figure  (i)  corresponds  to  the  east  leg,  (ii)  corresponds  to 
the  south  leg,  (lii)  corresponds  to  the  west  leg  and  (Iv) 
represents  the  north  leg.  Temperature  data  between  the  surface 
and  8m  depth  is  a  linear  interpolation  between  the  sea  surface 
temperature  and  the  first  CTD  Tow-Yo  data  bin.  Each  panel 
represents  the  observed  values  along  one  leg.  Each  panel 
starts  at  the  beginning  of  the  correspondent  leg  in  the 
following  directions  (NE-SE-SW-NW-NE) .  See  the  text  for 
further  interpretations. 
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Flgur*  I1I-C.3 


This  figure  shows  8  meter  vertical  shear  magnitudes  at  4 
equally  spaced  points  selected  along  the  west  leg.  The  strong 
shear  above  30m  depth  after  hour  150  in  the  profiles  of  point 
'600m'  is  simultaneous  with  the  rapid  variation  of  the  surface 
current  of  Figure  III-B.3.a,  and  with  the  attenuation  of  the 
rectification  process  at  the  thermocline,  as  shown  in  Figure 
III-A.5  on  the  profiles  at  point  '600m'  and  in  Figure  III- 
C.4.b  along  the  west  leg.  This  strong  shear  could  be 
associated  with  the  start  of  the  internal  mixed  layer  observed 
in  Figure  Ill-C.l.a  at  the  west  leg. 
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Profile  at  600  m 


Figure  III-C.3 
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Figure  IZI-D.l 


Current  meter  time  series  (a),  simultar'  ^th  the  MBAY 
survey  from  a  mooring  instrumented  with  temperature  (i), 
eastward  (U)  velocity  component  (ii)  and  northward  (V) 
velocity  component  (iii)  measurements.  The  hours  displayed 
along  the  x-axis  correspond  to  hours  after  0:00  of  September 
1,  1991.  The  MBAY  cruise  was  from  hours  120  to  186.  Spectral 
estimates  were  computed  for  128  hours  segments  and  are  shown 
in  Figures  (b) ,  (c)  and  (d) ,  for  temperature  (in  potential 
energy  units) ,  and  U  and  V  velocities  respectively.  The 
location  of  the  mooring  was  within  2000m  north  of  the  center 
of  the  noirth  leg,  as  shown  in  Figure  II-A.l.  A  semidiurnal 
tidal  signal  is  evident  on  the  temperature  and  U  time  series, 
which  is  nearly  aligned  with  the  canyon  axis,  though  the 
spectral  estimates  also  show  a  significant  peak  at  the  24  hour 
period.  The  temperature  time  series  show  significant  drops  of 
0.6^C  24  hours  apart,  though  smaller  ones  can  be  observed 
between  them.  The  persistence  of  the  spectral  estimates  in 
time  suggests  that  there  were  no  abnormal  features  present,  so 
the  present  data  can  be  assumed  to  be  typical  for  this  time  of 
the  year  in  this  area. 
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140  150 
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Potential  Energy  Power  Spectrum  for  128  Hr  subsets  (units  in  (m/s)^2) 
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Frequency 


U  Velocity  Power  Spectrum  for  128  Hr  subsets  (units  in  (m/s)^2 
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Frequency 


V  Velocity  Power  Spectrum  for  128  Hr  subsets 


Frequency 


rigure  ZXX-D-2 


Schematic  representation  of  vorticity  generation  mechanisms  in 
a  tidal  flow  over  topography: (a)  potential  vorticity 
conservation-column  stretching  and  squeezing;  (b)  differential 
friction  due  to  velocity  shear  (for  quadratic  friction  law) ; 
(c)  differential  friction  due  to  lateral  depth  gradients;  (d) 
asymmetrical  depth  profile  changes  [from  Robinson,  1985, 


Figure  III-D.2 
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Figure  IV-A.l 


Tidal  velocity  amplitudes  over  the  simplified  continental 
slope  taken  as  representative  of  the  topography  in  the  HBA¥ 
vicinity.  Figure  (a)  represents  the  cross-shore  (U)  horizontal 
tidal  component  increasing  in  magnitude  as  depth  shallows. 
Figure  (b)  shows  the  consequent,  bottom  induced  vertical 
velocity  (W)  amplitudes,  that  are  in-phase  with  the  horizontal 
cross-shore  velocities.  The  horizontal  velocity  amplitudes 
were  computed  assuming  kinetic  energy  conservation  with  a 
rigid  lid  asstimption.  The  vertical  velocities  were  then 
derived  by  continuity  assuming  no  velocity  normal  to  the 
vertical  boundaries.  This  velocity  field  was  used  for  the 
characterization  of  the  internal  tidal  field  in  the  M6AY  area 
and  for  the  study  of  nonlinear  internal  tide  propagation, 
where  it  was  taken  as  the  non-stationary  mean  velocity  field 
on  which  the  waves  propagate.  Figure  (c)  displays  the 
resultant  internal  tide  forcing  term  (N^/t>x)  as  derived  by 
Baines  (1982)  taking  into  accoxint  the  assumed  stratification 
for  the  MBAY  conditions,  as  discussed  in  the  text. 
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Figure  IV-A.l  (b) 
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Figure  IV-c.l 


Ray  paths  for  run  #1  along  the  grid  representing  an  attenuated 
continental  slope  relative  to  HBAY  conditions,  extending  for 
10  Km  between  the  18  Km  and  8  Kn  distances  with  depth  ranging 
from  200  to  1000m.  Figure  (a)  shows  vertical  sections  of  the 
rays  and  Figure  (b)  the  top  view  (horizontal  paths) .  See  text 
for  interpretation  of  these  results. 
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RUN  #1  -  RAY  PATHS  EVERY  2  Hr  OVER  A  TIDAL  CYCLE 


Figure  ZV-C.l  (b) 
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Figure  IV-C.2 


Ray  initial  velocity  amplitudes  for  run  #1  along  the  cross¬ 
shore  axis.  The  shelf  break  is  at  18  Km  and  the  rays  are 
propagating  towards  decreasing  values.  The  horizontal  axis 
identifies  each  ray  path  by  the  corresponding  starting  time. 
The  letters  'up*  and  'dw*  identify  initially  upward  and 
downward  rays  respectively.  The  velocity  amplitudes  were 
computed  from  the  wave  action  and  intrinsic  frequency  along 
the  rays  for  an  initial  wave  velocity  amplitude  of  0.15  m/s. 
This  run  was  forced  with  waves  having  six  times  the  tidal 
frequency  and  the  typical  length  scales  derived  in  IV-A.  See 
text  for  interpretation  of  these  results. 
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RUN  #1  -  RAY  Initial  Ve 


Figure  1V-C.3 


Ray  paths  for  run  #2  along  the  grid  representing  a  continental 
slope  similar  to  MBAY  conditions,  extending  for  5  Km  between 
9  Km  and  4  Km  distances  and  with  depth  ranging  from  200  to 
1000m.  Figure  (a)  shows  vertical  sections  of  the  rays  and 
Figure  (b)  the  top  view  (horizontal  paths) .  This  run  was 
forced  with  waves  having  twice  the  tidal  frequency  and  a 
typical  length  scale  of  7Km,  as  derived  in  IV-A.  See  text  for 
interpretation  of  these  results. 
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RUN  #2  -  RAY  PATHS  EVERY  2  Hr  OVER  A  TIDAL  CYCLE 


Figure  IV-C.3  (a) 
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Figure  IV-C.4 


Ray  initial  velocity  amplitudes  for  run  #2  along  the  cross¬ 
shore  axis.  The  shelf  break  is  at  9  Km  and  the  rays  are 
propagating  towards  decreasing  values.  The  x-axis  identifies 
each  ray  path  by  the  correspondent  starting  time.  The  letters 
'dw*  and  'up'  identify  initially  upward  and  downward  rays.  The 
eunplitudes  were  computed  from  the  wave  action  and  intrinsic 
frequency  along  the  rays  for  an  initial  wave  velocity 
eunplitude  of  0.15  m/s.  See  text  for  interpretation  of  these 
results . 
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RUN  #  2  -  RAY  Initial  Velocity  Amplitude  in  m/s 


Figure  IV-C.4 
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Figur*  IV>C.5 


Ray  paths  for  run  #3  along  the  grid  representing  a  continental 
slope  similar  to  NBAY  conditions,  extending  for  5  Km  between 
9  Km  and  4  Km  distances  and  with  depth  ranging  from  200  to 
1000m.  Figure  (a)  shows  vertical  sections  of  the  rays  and 
Figure  (b)  the  top  view  (horizontal  paths) .  This  run  was 
forced  with  waves  having  semidiurnal  tidal  frequency  and  a 
length  scale  of  7Km,  as  derived  in  IV-A.  See  text  for 
interpretation  of  these  results. 
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RUN  #3  -  RAY  PATHS  EVERY  2  Hr  OVER  A  TIDAL  CYCLE 


Figure  ZV>C.5  (b) 
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Figur«  IV-C.6 


Time  series  of  frequency  at  the  source  for  run  #3.  Figure  (i) 
corresponds  to  the  Eulerian  frequency  and  Figure  (ii)  to  the 
intrinsic  frequency.  The  horizontal  axis  represents  the  time 
during  a  typical  tidal  cycle.  The  effects  of  doppler  shifting 
are  evident  in  the  Eulerian  frequency  producing  a  modulation 
with  the  tidal  period;  however,  the  straining  effects  are 
visible  in  the  intrinsic  frequency  in  the  form  of  a  modulation 
by  a  frequency  twice  the  tidal,  canceling  the  symmetry  of  the 
Eulerian  frequency  modulation. 
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Figura  IV-C.7 


Evolution  of  frequency  along  each  ray  path  for  run  #3.  Figure 
(a)  corresponds  to  the  intrinsic  frequency  and  Figure  (b)  to 
the  Eulerian  frequency.  The  horizontal  axis  represents  the 
propagation  time  for  each  ray  since  it  started  from  the  source 
point  at  the  shelf  break.  Rays  are  marked  in  accordance  to  the 
time  on  the  tidal  cycle  at  their  beginning.  One  should  note 
the  range  of  the  intrinsic  frequency  extends  from 
approximately  the  inertial  frequency  (O.SxlO'^s"’)  to  3  times 
the  inertial  frec[uency.  The  Eulerian  frequency  shows  that  some 
rays  experience  a  strong  advection,  overcoming  their  group 
velocities  during  at  least  part  of  the  tidal  cycle,  as  can  be 
seen  by  the  negative  values. 
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RUN  #3  -  INTRINSIC  FREQUENCY 
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RUN  #3  -  EULERIAN  FREQUENCY 


Time  after  ray  start  (Hr) 


rigur*  zv>c.8 


Evolution  of  wavenumber  along  each  ray  path  for  run  #3.  Figure 
(a)  corresponds  to  the  horizontal,  cross-slope  wavenumber 
component  and  Figure  (b)  to  vertical  wavenumber  component.  The 
horizontal  axis  represents  propagation  time  for  each  ray 
since  it  started  from  the  scarce  point  at  the  shelf  break. 
Rays  are  marked  in  accordance  to  the  time  on  the  beginning  of 
the  tidal  cycle.  One  should  note  the  modulation  of  the 
wavenumber  by  the  spatial  derivatives  of  the  mean  flow.  The 
rays  that  start  during  the  flood  are  advected  onto  the  shelf 
where  no  changes  occur;  however,  when  they  propagate  back  to 
the  shelf  break,  though  the  vertical  component  has  significant 
variability,  the  horizontal  cos^onents  do  not  change  like 
those  starting  at  the  ebb. 
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RUN  #3  -  VERTICAL  WAVENUMBER 


Time  after  ray  start  (Hr) 


Figure  ZV-C.9 


Scheaatic  of  ray  paths  and  velocity  field  estimation.  A 
typical  characteristic  trajectory  (j)  starting  at  time  'tO' 
and  'x»z»»0'  in  the  space  time  domain  is  displayed.  The 
vari2d>les  't/'  correspond  to  the  instants  the  ray  arrives  to 
each  location  in  the  x-z  plane.  The  time  variant  velocity 
function  ' V* ,  starting  at  the  origin,  is  propagated  along  the 
ray  changing  its  parameters  in  accordance  with  changes  in  the 
mean  flow  and  location  in  space.  The  instantaneous  velocity 
estimates  at  time  'T*  are  then  computed  by  projecting  the 
characteristic  path  onto  the  plane  x-z  that  intersects  the 
time  axis  at  the  corresponding  time.  The  velocity  phase  lags 
are  obtained  from  the  time  lags  'T-t/'  and  the  time  the  ray 
crossed  the  origin. 
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Figar*  xy-c.9 
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Figure  IV>C.10 


Position  of  the  interface  of  a  two  layer  model  at  the  various 
phases  of  the  barotropic  tide  over  the  edge  of  a  slope.  Arrows 
denote  current  direction  [from  Baines,  1982  fig. 6]. 
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Figure  iv-c.ii 


Contours  showing  the  24  hour  velocity  averages  resulting  from 
model  run  #3  over  the  simplified  continental  slope.  Figure  (a) 
shows  the  onshore  velocity  component  (U) ,  and  Figure  (b)  shows 
the  alongslope  component  (V) ,  with  the  positive  direction 
being  into  the  page.  The  model  was  forced  with  semidiurnal 
internal  waves  at  a  point  source  at  25m  depth  in  the  shelf 
break.  The  parameters  for  this  run  are  summarized  in  Table  IV- 
C.l.  Strong  spatially  inhomogeneous  velocity  residuals  are 
observed  in  depth  ranges  from  30  to  70m,  though  smaller 
residuals  can  be  observed  in  the  profiles  at  7500m  distance 
between  150  and  200m  depths. 
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Figure  IV-C.ll  (a) 
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RUN  #3  -  24  Hr  Averaged  V  Velocity  In  m/s 


Figure  IV-C.ll  (b) 
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Figure  iv-c.12 


Con'tours  showing  the  velocity  profile  time  series  of  model  run 
#3  at  6500m  (a) ,  7500m  (b) ,  and  8250m  (c)  .  The  onshore 
velocity  component  (U)  is  displayed  in  Figures  (i)  and  the 
alongslope  component  (V) ,  with  the  positive  direction  being 
into  the  page  is  displayed  in  Figures  (ii)  .  The  model  was 
forced  with  semidiurnal  internal  waves  at  a  point  source  at 
25m  depth  in  the  shelf  break.  The  parameters  for  this  run  are 
summarized  in  Table  IV-c.l.  Wave-like  patterns  at  several 
frequencies  are  seen  in  all  figures,  though  the  profiles  at 
8750m  show  higher  frequencies  compared  to  the  others.  One 
should  note  the  24  hour  periodicity  in  the  V  component  at 
6500m,  between  40m  and  100m  depths,  that  can  also  be  seen  at 
50m  depth  in  the  profiles  at  7500m.  The  24  hour  averages 
displayed  in  Figure  IV-C.ll  show  strong  residual  currents  at 
these  locations  and  depths. 
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RUN  #3  -  U  Velocity  at  6500  m  (In  m/s) 


Figure  iv-c.12  (a) 
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RUN  #3  -  V  Velocity  at  6500 
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RUN  #3  -  U  Velocity  at  7500  m  (in  m/s) 


Figure  iv-c.12  (b) 
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Figure  IV-C.12  (c) 
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Figure  lV-C.i3 


Contours  showing  counts  of  hourly  occurrences  of  Richardson 
nuiBber  below  the  critical  value  of  0.25  over  a  24  hour  period. 
The  assumed  stratification  is  described  in  TeO^le  IV-c.l.  The 
high  number  of  occurrences  of  below  critical  conditions  at 
depth  ranges  of  50  and  100m  depth,  and  within  3000m  from  the 
shelf  break,  suggests  that  internal  tide  shear  instabilities 
might  exist  forcing  internal  mixed  layers  as  observed  during 
the  MBAY  experiment.  One  should  note  the  poor  resolution  of 
the  estimates  due  to  the  broad  resolution  of  the  ray  tracing 
numerical  method  that  strongly  enhance  the  edges  of  the 
internal  wave  beam.  Nevertheless,  the  relative  potential 
contribution  of  the  three  classes  of  rays  (deeper  rays 
starting  during  the  flood,  shallower  starting  during  the  ebb, 
and  intermediate  starting  during  slack  waters)  is  evident  in 
the  figure,  with  the  shallow  rays  displaying  a  larger  capacity 
to  generate  shear  instabilities. 


205 


HUN  #3  Hourly  Occurrences  of  \~\\  <  \)2.b  over  Z4  Hr 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

LH 

0 

LO 

0 

ir> 

0 

LD 

0 

LD 

! 

r— 

T— 

(NJ 

CVJ 

CO 

CO 

I  I  I  I  I  I  I  I 


J313UJ  Ui  L-ldaQ 


o 

a 

LH 

i 


Figure  IV-C . 13 
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